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(57) Abstract: A method creates an adjusted credit risk transition matrix and related matrices. The method uses the steps of obtain- 
ing a real probability measure (P) transition matrix, and adjusting the (P) transition matrix to be consistent with default probabilities 
under maitingale measure (0 by making the last column (default probabilities) in the (P) transition matrix consistent witii default 
^ probabilities calculated under (0). The method further scales the othw entries in the transition matrix (F) to compensate for ttie 
adjustment while retaining the relative weights among non-default classes firom die (JP) transition matrix in the adjusted transition 
matrix. A credit risk model for the analysis of complex credit instrument implements die method. The adjusted nunsiticHi matrix 
Tf[ may be for use in analysing complex credit instruments. The model may be implemented in computer software on a computer read- 
able medium. The adjusted transition matrix may be used to price credit swaps or value credit instruments. The model may have a 
^ simulation engine to generate valuation scenarios based on die adjusted transition matrix. The simulation engine may generate sce- 
narios based on a Monte Carlo mediod. The model may be extended for the risk analysis of complex portfolios of credit instruments 
Q with stochastic interest and foreign exchange rates. The model may have a Monte Carlo simulation engine to generate valuation 
scenarios based on the joint distribution of interest rates, spreads, foreign exchange rates and obligors* credit migrations consistent 
widi obtained or adjusted (/*) transition matrix. 
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CREDIT RISK ESTIMATION SYSTEM AND METHOD 
FIELD OF THE INVENTION 

The field of the invention relates generally to systems and methods for estimating credit risk 
associated with credit instruments and portfolios of credit instruments. 

BACKGROUND OF THE INVENTION 

Credit risk models are used to estimate the risk associated with credit instruments and 
portfolios of credit instruments, such as bonds and loans. Various models have evolved 
overtime. Many run simulations to generate possible valuation scenarios. Simulations are 
typically generated using a Monte Carlo method or variant. 

Credit risk models typically generate scenarios using probability distributions based on 
historical data. These probability distributions may be obtained in the form of a transition 
matrix from various sources, for example Standard & Poors™ and Moody's™. Such 
probability distributions are often referred to as {F), a real measure probability distribution. 

In credit risk models, transition matrices under (P) probability measure, or {F) transition 
matrices, are used to predict transition, or migration, from one credit class to another, "credit 
migration", over a given time period. A set of probabilities of migration, or transition, for 
different credit classes forms a transition matrix. 

Models are typically implemented in computer software as part of a credit risk system using 
coinpatible hardware. 
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Alternative matrices, models, systt ms and methods for estimating credit risk and credit 
migration associated with credit instruments are desirable. 



SUMMARY OF THE INVENTION 



In a first aspect the invention provides a method of creating an adjusted (P) transition matrix 
of probabilities of credit migration for different credit classes. The method uses the steps of 

• obtaining a transition matrix under probability measure (P) 

• adjusting the (P) transition matrix to be consistent with default probabilities under 
maningale measure (O) computed, for instance, from bond market prices by making a 
column of default probabilities in the (/') transition matrix consistent with default 
probabilities calculated under {(J) 

• scaling the other entries in the (P) transition matrix to compensate for the adjustment 
while retaining the relative weights among non-default classes from the (P) transition, 
matrix in the adjusted transition matrix. 

The method may adjust the (P) transition matrix in accordance with the following: 

R = {/^, Rf. } - are all credit classes 

Jif, - corresponds to default state 

T'' = : j(to*f)j - adjusted (P) transition matrix for the time interval [/o,^] 
P,.j(^o*^\J - transition probability from credit class R. to credit class Rj 

p-'^'^-'^'' \v,p:,. j<K f"^^^ 

PkA'o>U) = OiT J<K and p^.,: (/o,/.) = 1 
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1 - A.A- 

A.A' - obtained (P) probability of default for the credit class 

(from the obtained (P) transition matrix) 

- (Q) probability of default for the credit class R 

G - a set of credit classes (can be the whole set R = {/i; R^] 

or some subset 

I ^ ' ' denotes the number of elements in G , if G is empty, put 

c; = o 

N{) 

- standard normal cumulative density function. 

Alternatively, the method may adjust the (/>) transition matrix in accordance with the 
following: 

C = II A?, . ) II - adjusted (/') transition matrix for the time interval [/„,/,] for 
obligor n 

P".t('oA) - obligor's /; transition probability to migrate from credit class 

R, to credit class 

PK.jito,'i) = 0 \f j<K and p;.^-(/o./,) = l 
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^0 ^0('o»'l>'~ 1 7 Y 



- obtained (P) probability of default for the obligor n over 
interval [/q,/, ] if it is in the credit class R, at time /„ 

</" - (0 probability of default for the obligor /; over interval 

if it is in the credit class R at time 

r, - instantaneous risk-free interest rate at time / 

y," - obligor n stock log-return over interval [/q,/,] 

a," - obligor // drift function of time f 



In a second aspect the invention provides a credit risk model for the risk analysis of complex 
portfolios of credit instruments. The adjusted (P) transition matrix may be used to value 
credit swaps and other complex credit instruments. The model may have a simulation engine 
to generate valuation scenarios based on the adjusted (P) transition matrix. The simulation 
engine may generate scenarios based on a Monte Carlo method. The model may be 
implemented in computer software on a computer readable medium. 

In a third aspect the invention provides a method of pricing a credit swap. The method 
utilizes a number of steps: 
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• {P) transition matrix over a risk horizon [/„./, ] and the (0 default probabilities over the 
same period [/„,/,] are obtained 

• an adjusted {P) transition matrix and (0 transition matrix over /, ] are created 

• a first series of forward (0 transition matrices for standard observation times 
/ = /|./,,A3,... IS built from the (0 transition matrix and an initial term structure of 
default probabilities under O over [/„,/,] 

• the first series is translated to produce a second series, having premium payment dates 
as the observation times by working with 'fractional' powers of transition matrices for 
subintervals 

• a third series of transition matrices is derived from the second series for co-migration on 
credit classes for pairs of obligors 

• a risk-free discount factor and probabilities contained in the third series are used to 
discount cash flow for each payment date over the period, and the discounted cash 
flows are summed for all payment dates. The sum is subtracted from a payoff function 
for fixed payments to give the credit swap price. 

In a founh aspect the invemion provides an extended credit risk model for the risk analysis of 
complex portfolios of credit instruments with stochastic interest and foreign exchange rates. 
The model may have a Monte Cario simulation engine to generate valuation scenarios based 
on the obtained or adjusted {P) transition matrix. The simulation engine may generate 
scenarios from a joint lognormal distribution of the following risk factors: 

• country/industry market indices 

• obligors' idiosyncratic components 
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• base forward rates /Um/i)./(/,./:),...../V/,-i./J for time intervals [/o./J, [/,./J [^..p/J 

corresponding to given maturities /„./, t^, 

• credit spread .v,(/,) for the highest credit class and for maturity equal to the first 
bucket point f, on the corporate zero curve for each currency 

• incremental spreads A.v (/()= .v {/,)-. v ,(/,), / = 2 A' - I , for maturity A (by credit class 

and currency) 

• the difference Av^. ^(/J=.v^. ^(/J-.v^. ^(/j) between the spreads for the lowest credit 
rating for the longest and shonest maturity (by currency) 

• foreign exchange rates. 

From risk factors logrcturns for each simulation run the model may compute: 

• obligors' credit ratings at risk horizon based on country/industry market index logreturns 
and idiosyncratic components 

• base zero interest rates for all currencies based on forward rates 

• credit spreads .v (/j ) = v ,(/,) + A.y, (/, ) for maturity /, . / = 2,. ..A' - 1 (by currency) 

• credit spread .v^. /i^j^s^. ; + A.v^. ^(/^) for the lowest credit rating and longest 
maturity (by currency) 

• the rest of the spreads for maturity by partitioning the interval [O.A;. 

proponionallv to .v, (/, ) bv taking a f/^ ) = sjt^ jl^Ljihl (by currency) 

• spreads for other maturities , by linear interpolation between .v. (/, ) and s^{t^) 

(by currency) 
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• corporate zero curves at risk horizon by aoding the appropriate spread to the base curves 
(by currency) 

• foreign exchange rates. 

In other aspects the invention provides various models, matrices, software, systems and 
methods incorporating combinations and subsets of the aspects described above, or utilizing 
the aspects described above. 

BRIEF DESCRIPTION OF THE DRAWINGS 

For a better understanding of the present invcmion and to show more clearly how it may be 
carried into effect, reference will now be made, by way of example, to the accompanying 
drawings which show the preferred embodiment of the present invemion and in which: 

Fig. 1 is a block diagram of the structure of a credit risk model according to 

the preferred embodiment of the invention. 
Fig. 2 is an example of Portfolio Summary (analytic engine was not used) 

output from the credit risk model of Figure I. 
Fig. 3 is an example of Obligor Summary output from the credit risk model 
of Figure 1, 

Fig. 4 is an example of Distribution of the Portfolio Forward Value output 

from the credit risk model of Figure 1, and 
Fig. 5 is an example of Marginal Risk vs. Exposure output from the credit 

risk model of Figure 1, and 
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Fig. 6 is a block diagram of the structure of a credit risk model with 
stochastic interest and foreign exchange rates according to the 
preferred embodiment of the invention. 

DETAILED DESCRIPTION OF THE PREFERRED EMBODIMENTS 

Probability distribution (probability measure) (P) is a probability distribution of risk factors 
based on historical data. In credit risk models, probability distribution (P) is used to predict 
transition, or migration, from one credit class to another, "credit migration" A group of 
probabilities of migration, or transition, for different credit classes forms a transition matrix. 
A series of transition matrices can be built for different time intervals and for different 
obligors. 

In the preferred embodiment, a given, or obtained, transition matrix under the probability 
measure (P) is adjusted to be consistent with default probabilities under martingale 
probability measure (Q). In the preferred embodiment, the (P) transition matrix is adjusted 
by making a column of default probabilities consistent with default probabilities calculated 
under (0. The other entries in the transition matrix are scaled to compensate for the 
adjustment while retaining the relative weights among the non-default classes in the transition 
matrix. 

The adjusted (P) transition matrix may be incorporated into a credit risk model to analyze 
complex portfolios of credit instruments, including credit derivatives. For example, the 
model may use the adjusted (P) transition matrix in place of the given (P) transition matrix in 
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a simulation engine to generate scenarios. The simulation engine may generate scenarios 
based on any Monte Carlo method. 



The model may be implemented in computer software as part of a credit risk system using 
compatible hardware. 



The adjusted (P) transition matrix may be used in the pricing of credit swaps. Credit swap 
pricing may be incorporated into the credit risk model as it is consistent with the 
methodology used in the model. 

This description has four Parts. 

Part 1 describes assignee s basic credit risk model. CreditVaR'-^'. which utilizes deterministic 
interest and foreign exchange rates. Each obligor s standardized equity returns are 
decomposed into weighted average of market indices returns (multi-beta model) where the 
weights are specified to appropriately reflect the obligor s participation in the corresponding 
markets and to model obligor s idiosyncratic returns. Part I has Appendix I, 2 and 3. 

Part 2 describes the derivation of an adjusted (P) transition matrix and related matrices. In 
particular the relationship between P and O probabilities for obligors to migrate from one 
credit class to another including probabilities of default is derived. An adjusted transition 
matrix and related matrices are created based on the derived relationship. The adjusted (P) 
transition matrix and related matrices allow portfolios of different instruments, including 
credit derivatives, to be analyzed. In the preferred .mboc;imem, the adjusted (P) transition 
matrix and related matrices have been incorporated into an extended version of the credit risk 
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model of Part I. Part 2 also describes a method of pricing credit swaps using the adjusted {P) 
transition matrix and related matrices. In the preferred embodiment, the method has been 
incorporated into the credit risk model of Part I. The adjusted transition {P) matrix, related 
matrices and correlations between equity returns are used to compute joint probability 
distribution of obligors' credit migrations. 

Part 3 describes an extended version of the credit risk model that allows for stochastic interest 
and foreign exchange rates. The model generates correlated interest rates, foreign exchange 
rates and obligors' credit migrations according to their joint probability distribution. As 
stochastic rates are permitted valuation must be performed for each instrument in a generated 
scenario. 

It will be recognized that the adjusted {P) transition matrix and related matrices, pricing 
method and stochastic features may also be incorporated into other compatible credit risk 
models. - 
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Part 1 - Basic Credit VaR Methidoio^y 



1. Introduction 



CreditVaR'^' is assignee s proprietary model for measuring and analyzing credit risk in a 
portfolio context. This description is made in reference to Credit VaR as the preferred 
embodiment of the credit risk model. Those skilled in the art will recognise that many of the 
principles described herein are equally applicable to other credit risk models, such as J. P. 
Morgan's implementation ofCreditMetrics™. 



The methodology used in CreditVaR is based on credit migration analysis and on the 
assumption that an obligor s credit migration is driven by its asset value. For additional 
discussion of credit migration analysis, please see CreditMetrics'^' - Technical Document, 
April 2, 1997, J P. Morgan & Co. Incorporated that is hereby incorporated by reference and 
available at www.creditmetrics.com. 



The methodology structure and implementation of CreditVaR are shown in Figure 1. 
Information required for credit risk analysis includes. 



Obligor data: 

Information about obligors is organized into a database containing details of their 
credit ratings, industries, and countries. 



Portfolio data: 
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Information about financial positions is organized into portfolios of exposures. It 
allows to cover different types of instruments such as fixed income instruments, loans, 
commitments, letters of credit, etc. 

Market dcitci^ transition probabilities and correlations: 

These include yield curves, spread curves, foreign exchange rates, transition 
probabilities from one credit rating to another for different credit rating systems 
(Moody 8 states. Moody 18 states, S&P 8 states), correlations between market 
indices 

The CrcdiiVaR methodology may be implemented both analytically and as a Monte Carlo 
simulation. It calculates two risk measures: standard deviation (in the analytic and simulation 
engines) and percentile level (in the simulation engine) of the portfolio value distribution for 
a given time horizon. The regulatory capita! calculation is based on the percentile level, 
therefore, only simulation engine is used for the capita! calculation. Accordingly, in the 
preferred embodiment the analytic engine was removed. Following sections describe the 
methodology and a possible implementation in more details. 

2. Risk Slatistics and IVIeasurement 

The CreditVaR methodology can assess the impact of changes in debt value for a given 
portfolio due to credit quality movements of issuers (obligors) - including downgrades, 
upgrades and possible defaults - which can occur within the time horizon, typically a period 
of one year. 
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There are two computational engines developed in CreditVaR: Analytic Engine (AE) and 
Simulation Engine (SE). Both are capable of producing mean and standard deviation for the 
portfolio's forward value (risk horizon later). The two sets of numbers should be quite 
comparable as the same modelling framework is embodied in both engines. Though 
distribution of value for a portfolio of bonds is typically very much skewed, these statistics 
contribute to some understanding of credit risk present in the portfolio. 

The Simulation Engine generates a distribution of forward portfolio values, and therefore has 
the advantage of estimating any statistic, including ones deemed useful in risk measurement. 
In panicuiar, i 'aluc-cu-Risk for any percentile level, together with a confidence band, can be 
computed. (Another risk statistic that can be produced from SE is average shortfall, which 
is the expected loss given that losses exceed a given level The expected cxcession of a 
percentile level, which is the expected loss given that the loss is more extreme than the given 
percentile level, is another very useful risk measurement that can be obtained using SE 
engine.) Results become more accurate as number of simulations increases, but will exhibit 
cost in terms of time and memory as with any Monte Carlo method. 

With additional computation, both engines AE and SE produce marginal standard deviation 
by each obligor, i.e., the difference between the standard deviation of the entire portfolio and 
that of the portfolio ex the obligor. Plotting marginal standard deviation expressed as 
percentage of mean value for the given asset against the exposure size for each obligor, it can 
be used to identify that pan of the ponfolio which has concentration of risk, large exposure 
size and high percentage of marginal standard deviation. All of these measurements are 
useful in risk monitoring, as well as for capital allocation purpose. 
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The current CreditVaR method can be extended to produce measurement of risk which are 
due to obligor-specific returns. Risk due to market indices, in fact to any subset of market 
indices, can also be obtained. 

3. Pre-processing 

The methodology requires that the three types of data mentioned in introduction be given as 
inputs. They are then used to prepare the portfolio in an input format applicable to the 
CreditVaR Enuine(s). We now indicate the required pre-processing. 

For each instrument in the portfolio, there should be sufficient information to calculate its 
forward value for any credit migration with respect to a credit rating system chosen for 
CreditVaR. Specifications such as coupon rate, term to maturity, yield curve, spreads, etc. 
are needed. Recovery rate and its standard deviation, i c. the seniority of the debt, are also 
required in case default occurs. We assume that obligor's exposures can be aggregated or 
netted across instruments having same seniority. 

The credit quality movements can be tabulated in a transition matrix (for the given risk 
horizon) of probabilities for an obligor in one rating to end up in other ratings, including 
default, or remain unchanged in its rating An additional ingredient essential in the pre- 
processing is a module which quantitatively measures the portfolio effect of credit by 
accounting for correlations of asset returns for the obligors. Once equipped with these asset 
correlations (which are on average typically between 20% and 35%), joint rating changes can 
be modelled across these obligors. 
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Ultimately, it turns out that in t lis model the only parameters which affect the risk of a 
portfolio are the two already men.ioned. namely, the transition probabilities for each obligor 
and the correlations between assci returns. In this framework, the formulation of risk 
estimation is then reduced to considering the ^ standardized" (mean 0, standard deviation 1) 
asset returns, where the only parameters to be estimated are these correlations. 

It should be noted that interest rates and spreads are assumed to be deterministic in this Part 
1 In Parts 2 and 3 interest rates and foreign exchange rates are allowed to be stochastic. As 
stochastic rates are permitted valuation must be performed for each instrument in a generated 
scenario. Part 3 describes a method and means for valuing such scenarios. 

As asset values are in general difficult to observe, we use instead the correlation between 
equity returns as its proxy. The obvious advantage of this is that equity data are more readily 
available than data of credit spreads or of actual joint rating changes. In order to keep the 
size of correlation matrix reasonably small to be useftil in computation as well, the 
methodology relies on some form of component analysis. Basically, each obligor's 
standardised asset returns are decomposed into a weighted average of some 'benchmark' 
country-industry indices where the weights are specified to appropriately reflect the obligor's 
participation in the corresponding markets, plus a residual component that is associated to 
obligor's idiosyncratic returns. These weights are specified for each obligor as part of input 
stored in Obligor database. (An example of such specification would be: Company X 
participates 70% in German Chemicals and 30% in German Electronics, and 15% of the 
movements in X's equity returns are company-specific.) In this way, the methodology 
incorporates the mix in obligor's market capitalization as well as obligor-specific movements 
in returns when computing credit risk of the portfolio. Given these weights and a correlation 
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matrix for the indices, one can easily obtain a matrix of correlations between ail of the 
obligors present in the portfolio. Appendix 2 presents in detail this methodology. 

We remark that in applying CreditVaR to a number of ponfolios there arc several 
computational modules which can be made external to the pre-processing of individual 
portfolios. For example, the asset correlations between obligors can be stored in a database, 
and each time a ponfolio is encountered for credit risk computation selecting appropriate data 
can be considered as part of the pre-processing. 

4. Analytic Engine (AE) 

Mean for the portfolio's forward value can be calculated in a straightforward way, using 
transition probabilities and forward bond prices which have already been valued in the pre- 
process. Computing standard deviation requires more work, but it reduces to computing 
standard deviations of various obligor-pair subportfolios. The latter is straightforward as well 
since joint probabilities of credit quality co-movements arc made available in the pre-process. 
iVlarginal standard deviations can be computed in a similar fashion. 

5. Simuiation Engine (SE) 

First, threshold levels of standardised asset returns representing credit rating changes are 
determined for each obligor in the portfolio using the transition probabilities: this part 
depends only on the obligor's credit class and not on the obligor itself. Using the asset 
correlations as the covariance matrix, samples of standardised multivariate normal random 
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variables are generated. These are then bucketed using the threshold levels, thereby 
generating scenarios of credit rating co-movements for the obligors. 

For each of these scenarios, portfolio is revalued quickly using pre-processed prices of 
exposures. Each time a default occurs in scenario, i.e.. when the sampled asset return value 
of an obligor is below the default threshold level, a random recovery rate is generated 
according to a beta-distribution whose defining parameters are governed by seniority 
associated to the obligor s bonds Finally, we obtain a distribution of portfolio values and 
from it the relevant risk statistics. 

6. Implementation 

A CreditVaR software program utilizes the methodology described above. The core 
computational engine is implemented in C++'^'. The imerface for Windows''^' could be a 
combination of Microsoft Access'^^'/Microsoft Excel™ Windows version will consist of a 
Microsoft Access/Excel application and a DLL with the compiled computational engine. The 
selection of hardware is at the discretion of the user provided it has sufficient processing 
capability and is compatible with the programming environment chosen. Other 
programming environments could be used as would be evidem to a person skilled in the art. 

CreditVaR program uses a country-industry covariance matrix to compute asset correlations 
between the obligors in the portfolio Matrices may be comained on a computer readable 
medium, for example on a hard disk, floppy disk, optica! disk, or random access memory, as 
appropriate for the particular computing environment anc. stage of use. Matrices may be 
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contained in an appropriate signal for transmission, such as across the internet or another 
computer network 

Computation engine consists of 2 parts - analytic and simulation. Analytic solution produces 
total mean and total and marginal standard deviation of every obligor. Simulation produces 
total mean and standard deviation of every obligor, as well as Value-at-Risk numbers for the 
whole portfolio. For initial implementation time horizon is taken to be 1 year. The model can 
be expanded to handle other credit instruments like loans, loan commitments, etc. 

Inputs 

The current implementation of the Global Analytics CreditVaR model evaluates the credit 
risk of a bond portfolio based on the following information: 

1 . Obligor data 

• Current credit ratings of the issuers (obligors), company-specific risk of the obligors, 
capitalization data 

2. Market Data 

• Current zero-rate curves and current spread curves corresponding to different credit 
ratings and currencies 

• Forward ( I year from now) zero-rate curves and forward spread curves corresponding 
to different credit ratings and currencies 

• Transition probabilities matrix for credit rating migrations 

• Covariances between the market indices underlying obligors' capitalization 

3. Spot and forward (1 year) exchange rates 

4. Market Data exposure data (PortEfolio) 
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• For bonds: principal amou n of the bond, coupon rate or spread, seniority (defined by 
recovery rate and its standard de\'iation), maturity of the bond, currency, etc. 

• For loans and loan commitments: total authorized amount, loan amount, letters of 
credit amount, coupon or spread, fees on letters of credit and unused portions, 
seniority, etc. 

Outputs 

The program produces the following outputs: 
I Aggregated by obligor: 

• Total present value of the bonds by a given obligor 

• Total expected I year forward value of the bonds (obtained analytically and by 
simulation) 

• Total standard deviation of 1 year forward value of the bonds (obtained analytically 
and by simulation) 

• Stand-alone VaR at different confidence levels 

• Marginal standard deviation and marginal risk of a given obligor 

• Delta standard deviation and delta VaR 

2. Totals for the portfolio: 

• Total present value of the portfolio 

• Total expected portfolio value in I year from now (obtained analytically and by 
simulation) 

• Total standard deviation of the portfolio value in 1 year from now 

• Skewness and kunosis of the portfolio distribution 
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• Percentiles of the portfolio value (VaR numbers of the portfolio) and expected 
excession at 0. 1%, 0.5%, 1%, 5% and 10% levels 

• 95% confidence intervals for the VaR values 

Marginal standard deviation and marginal VaR for a given obligor / are defined as 

where S'ri\^^^^^ ^'<^'/Cu/ ' are portfolio standard deviation and VaR, and STD . 

' ^K.iu/ uhhyor, - are standard deviation and VaR of the portfolio without positions 

corresponding to this obligor. 

Obligor's Delta standard deviation (Delta VaR) is defined as the sensitivity of the ponfolio^s 
standard deviation (VaR number) with respect to the aggregated obligor* s / position value 
Pubu^or, > multiplied by : 

^Publtxior, 

DeltaSTD^,^^^^^ and Dc/tal^aR^^^^^^^ satisfy following equations [3]: 

In other words. De/taSTD^,,^,,^ (DeliaVali^,,^^ ) is an obligor / contribution into the 
portfolio's total standard deviation (VaR). 
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Appendix 1 to Part 1 

Sample Input and Output Data 

Obligor data table (example) - contains the following information: 

• Short and Full name of the obligor 

• Current credit rating of the obligor 

• Company specific % - the percent of the company's risk that is not explained by 
company's sensitivity to a particular country or sector of industry 

Country I (2, 3) and Industry 1 (2. 3) for each country and appropriate percentages - 
capitalization of the company 
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Short 
Name 


Long Name 


Rating 


Company- 
Specific % 


Countr 

y 1 


Capitalization 
in Country 1. % 


lndLj<>trv 
1 


in Industry 1 , 


AAAAKOT 


General Electric 


Aaa 


20 


us 


100 


GNRL 


100 


AAAABPD 


Hydro-Quebec 


Aa 


61 


CA 


100 


GNRL 


100 


AAAAAUD 


Ontario, 
Province Of 


Aa 


72 


CA 


100 


GNRL 


100 


28664 


Dharmala Sakti 
Sejahtera. PT 


Ba 


51 


ID 


100 


GNRL 


50 


29715 


Citicorp 


A 


53 


US 


100 


GNRL 


20 


29841 


Anwar Sierad, 
PT 


B 


51 


ID 


100 


GNRL 


100 


27508 


BOC Group Pic, 
The 


A 


52 


GB 


70 


GNRL 


50 


29842 


Dharmala Inti 
Utama 


Ba 


51 


ID 


100 


GNRL 


100 



Short 
Name 


Long 
Name 


Capital- 
ization 
in 

Industry 
1.% 


Indust 
ry 

2 


Capital* 
ization in 
Industry 
2,% 


Industry 
3 


Capital- 
ization 
in 

Industry 
3.% 


Count 
ry 
2 


Capital- 
ization 
in 

Country 

2.% 


AAAAKOT 


General 
Electric 


100 














AAAABPD 


Hydro- 
Quebec 


100 














AAAAAU 

D 


Ontario. 

Province 

Of 


100 














28664 


Dharmala 
SaKti 

Sejahtera, 
PT 


50 


AUTO 


20 


FOOD 


30 






29715 


Citicorp 


20 


BFIN 


80 










29841 


Anwar 
Sierad. PT 


100 


ELCS 


50 






US 


30 


27508 


BOC 

Group Pic. 
The 


50 














29842 


Dharmala 
Inti Utama 


100 
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Current zero-rate curves for compiting present values of the bonds ffragmem): 



Currency 



AUD 



AUD 



AUD 



AUD 



AUD 



AUD 



AUD 



AUD 



AUD 



BEF 



BEF 



BEF 



Maturity 



BEF 



BEF 



BEF 



BEF 



BEF 



BEF 



CAD 



CAD 



CAD 



CAD 



CAD 



CAD 



CAD 



CAD 



CAD 



10 



Aaa 



Aa 



Baa i 



Ba I 



Bl 



0. 053250 I C . 0 547.66 I 0.054966 I 0 056666 I 0.061166 I Q"Q fi7P^ 
' O.C53^34 I 0.054536 f 0.056133 i 0.062178 i Q nfi7ftfi7 



Caa 



0.054404 



0.057863 



0.062668 



0.077766 
0.081702 



! Q Q^^Q^^ > Q 058845 I 0 .065964 I 0.072 3;>4 i o nQc^non 



0.065801 



15 



20 



30 



10 



0.068452 



0.059883 I 0.06 0 592 i 0.062078 | 0.071337 I 0.07853.S 1 n 'moPAy 



0 062786 I 0.063617 I 0 065720 I 0.075229 



0.065865 I 0.066699 I 0.069597 j 0.079848 



0.070039 I 0.068845 I 0.069800 I 0^072824 I Q Qfid.-^r^:^ 



0 070246 I 0.077058 I 0.078239 I 0.081695 I Q OQ'^QQd 



0.084363 1 0.125788 



0.090017 



0.094863 



0.108289 



0.068616 1 0.075 739 1 0.077275 1 0 081163 1 0;10Q634 I n iic;7io 



0.044/UO I 0.040700 i 0.041400 1 0.043100 | 0.047600 I 0053700 



0.145838 



0.160516 



0.160516 



0.160516 



0.047977 \ 0.042845 I 0.043648 I 0.045248 I 0 051 292 I Q.'osSQQI 



0 051217j 0.045657 
0.055735 i 0.050484 



0.046252 I 0.048066 I 0.055175 I 0.061536 
0 0511951 0. 052688 | 0.061 934 I 0.069137 



15 



0.059462 I 0.055583 I 0.056418 I 0.058528 1 0.068088 I 0.077251 



0.063016 f 0 059309 I 0.060148 I 0.063046 I 0 0 73:^51 n nfi-^^^o 
Q nfifi97fl I n nf^oAif^ I n r>goo-r^ . ^ : u.uoooo^ 



30 



10 



15 



20 



30 



0 066228 I 0.062416 I 0 063372 I 0.066395 j 0.077919 I n 



0.064200 



0.070827 



0.084151 



0.100376 



0.118753 



0.139378 



0 068648 I 0.0651 0 -rra ^2U9 I 0 0 69446 I 0.0827Q 1 ) n nZn^^ ■ n o.n.. 
0 .062548 I 0.0 69235 I 0 07-0766To-Q7;65^ 0.094125 I a?09.".^ r^4^ 

0.0543081 0.058808 I 0 064QnR I 0 075408 

0.079768 



0.046084 1 0.051908 \ 0.Q526Q8 



0.0480691 0.051797 0.0526001 0.054197 1 0.060242 1 0.Q659.T1 



0.051818 I 0.05 3409 | 0 054000 ) 0.055812 ( 0 062919 I 0.069277 



0.058771 



0.062437 



0.067168 



0.069212 



0.070499 



0 055384 I 0 056091 I 0.Q57S74 



0.057811 [ 0.0586381 0.0607P7 



0 0601 02 I 0.060929 I 0.063788 



0.064606 



0.066522 



0.065561 i 0.068582 



0.0676191 0.070830 



0.066785 I 0.073954 



0.070182 



0.073943 



0.080094 



0.083975 



0.067655 I 0.069103 ( 0.072814 I Q nQORTT 



0.079251 



0.083994 



0.090640 



0.095206 



0.105090 



0.120275 



0.138775 



0.156078 



0.187307 



CHF 



CHF 



CHF 



CHF 



CHF 



CHF 



CHF 



CHF 
DEM 



0.029251 



0.029647 



0.029943 



0 018100 ) 0. 0188001 0.020500 



0.021334 I 0.022137 1 0.0 23738 
0.024 503 I 0.025101 f 0 



.026913 



0 030250 I 0 030960 ( 0 032465 



0.025000 



0.029770 



0.033995 



0.041641 



10 



0.0304721 0.035345 i 0.036174 I 0.038262 I 0 0477^^ 



0.031100 



0.035476 



0.040351 



0.048827 



0.041600 



0 049285 



0.062777 



0.079720 



15 



0.031O12 I 0. 040189 ! 0 041027 I 0.043866 I 0 n.S4 l ^R 



20 



0.031123 1 0.044627 1 0.045579 1 0.048570 



! 0.056809 I 0.097483 
I_a064209 I 0.118510 



30 



0 031259 [ 0 047376 I 0.048456 i n nqi^ 



0 031582) 



0.0600251 0.070555 I 0.1 34308 



0.055377 
0. 040200 f 



0.056900 

0.040900 t 0.0426QQ 



0.064593 



0.075804 



0.162180 



I 0 060904 I 0.080289 I 0.09627fi 1 n ir57^ 



DEM 



DEM 



DEM 



DEM 



DEM 



DEM 



DEM 



DEM 



0.043600 i 



0.045953 I 0 041731 I 0.042534 I 0 0441 33 
0050143 I 0.04527 9 I 0J)A^R7Alc^nAy^ 



0.054534 I 0.049986 I 



0.058249 I 0.054607 I 



10 



0.050697 I 0.052191 
0.0554 41 I 0.057548 



0.047100 



0.050175 



0.054801 



0.061435 



0.053200 



0.055873 



0.063700 



0.069703 



0.061166 I 0.083789 



0.068639 



0.067092 



15 



0 061 660 I 0T058466 I 0.059305 I 0.062197 1 0 nv^AOA 



0.065009 



20 



066847 



30 



0.062601 
0.064869 



0.063563 I 0.066605 I 0.078204 



0.065974 I 0.0692151 0.082483 



0.076240 



0.082689 



0.088840 



0.093849 



0.099872 



0.117637 



0.138390 



0.154962 



.070431 I 0.0/4380 I 0.076066 1 0.080.398 i Q.iojss; i 0 1^20891 i g 
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F'orward ( I year) zero-rate curves tor computing forward values ot the bonds f fragment): 



Currency 


Maturity Aaa 


Aal Al Baal Ba 1 Bi Caa 


AUD 


U 


0,055559 


0.053202 { 0.054107 1 0.055601 1 0.063192 t 0.068470 1 0 085654 


AUD 


1 


0.055559 


0.053202 1 0.054107 ! 0.055601 i 0.063192 1 0. 068470 1 0.085654 


AUD 




0.060177 


0.057528 1 0.058066 1 0.059937 | 0.068371 1 0.074863 I 0.103750 


AUD 


A 

4 


0.065036 


0.061292 1 0.062003 1 0.063435 10.073895)0.0813701 0.118040 


AUD 


6 


0.067908 1 0.064212 1 0.065066 1 0.067236 1 0.077591 1 0.087239 \ 0 133997 


AUD 


n 

y 


0.0701 55 


0.067162 1 0.068010 1 0.071044 1 0.081944 1 0.092575 1 0.153663 


AUD 


1 4 


0.071248 


0.069894 1 0.070867 1 0.073988 1 0.086007 j 0.096861 1 0.166664 


AUD 


\ 9 


0.071 1 48 


0.078271 1 0.079478 I 0.083029 1 0.097859 I 0.110491 1 0.165043 


AUD 




0.0691 50 


0.076488 1 0.078053 1 0.082017 1 0.1020201 0.1174221 0.163480 


BEF 


0 


0.051264 


0.044994 1 0.045901 | 0.047400 


1 0.054998 1 0.060293 1 0.077496 


BEF 


1 


0.051264 


0.044994 1 0.045901 1 0.047400 


1 0 054998 1 0.060293 { 0.077496 


BEF 


2 


0.054490 


0 048145 1 0.048686 1 0.050557 


1 0.058983 1 0.065477 | 0.094266 


BEF 


4 


0.05851 1 


0.052945 1 0.053657 1 0 055099 


10 065547 10.073032 1 0.109610 


BEF 


6 


0.061942 


0.058084 1 0.058942 1 0.061 122 


10.071542 10.081230 1 0.128113 


BEF 


9 


0.065071 


0.0613971 0.062252 1 0 065285 


0 076252 1 0.086928 1 0.148053 


BEF 


14 


0.067782 


0.063984 1 0.064959 1 0 068079 


) 0.080118 1 0.091007 1 0.160574 


BEF 


1 9 


0.069923 


0.066404 1 0.067530 1 0.070851 


10.08458110.095219 1 0.194843 


BEF 


29 


0.063169 


0.0702331 0.0717931 0.075803 1 0.095765 1 0.1115061 0.192460 


CAD 


0 


0.050057 


0.051687 


0.052592 1 0.054087 I 0.061679 1 0.066959 I 0.084146 


CAD 


1 


0.050057 


0.051687 


0.052592 1 0.054087 j 0.061679 j 0.066959 1 0.084146 


CAD 


2 


0.054696 


0.054160 


0.054697 1 0.0565661 0.064981 1 0.071461 I 0.100264 


CAD 


4 


0.057776 


0.056255 


0.056963 1 0.058392 


0 068789 10.076227 1 0.112638 


CAD 


6 


0.060901 


0.058799 


0.059647 t 0.061801 


0.072090 10.081660 1 0.127933 


CAD 


9 


0.064270 


0.061017 


0 061858 1 0.064847 


0.075638 10.086136 1 0.146042 


CAD 


1 4 


0.068690 


0.065519 


0.066492 1 0.069609 


0.081631 1 0.092501 j 0.162067 


CAD 


1 9 


0.070443 


0.067297 


0.068415 1 0.071707 


0.085316 1 0.096824 1 0.193509 


CAD 


29 


0.071351 


0 068202 


0.069677 1 0.073458 


0.09200010.1062511 0.191367 


CHF 


0 


U .U^4lO I D 


0.024579 


0.025484 j 0.026986 


0.034562 1 0.039871 | 0.057028 


CHF 


1 


0.032815 


0.024579 1 


0.025484 1 0.026986 


0 034562 1 0.039871 1 0.057028 


CHF 


2 


0.031626 


0.027719 } 


0.028266 1 0.030135 


0.038522 10.045007 1 0.073527 


CHF 


4 


0.031006 


0.033310 1 


0.034023 i 0.035479 


0.045844 1 0.053306 1 0.089466 


CHF 


6 


0.031270 


0.038247 { 


0.039099 1 0 041253 1 


0 051598 10.0611551 0.107084 


CHF 


g 


0.031604 


0.042673 1 


0.043527 1 0.046495 1 


0 05740110.067953 1 0.127399 


CHF 


14 


0.031511 


0.046548 1 


0.047518 1 0.050605 j 


0.06257210.073430 1 0.141238 


CHF 


19 


0.031553 


0.048939 1 


0.050041 1 0.053304 j 


0.066718 10.078210} 0.168899 


CHF 


29 


0.031785 


0.056686 1 


0 058239 1 0 062326 1 


0.082248 1 0.098595 i 0.166578 


DEM 


0 


0.048311 


0.043265 1 0.044171 I 0.045669 1 0.053260 I 0 058552 1 0 075739 


DEM 


1 


0.048311 


0.0432651 0.044171 1 0.045669 I 0 053260 1 0.058552 1 0 075739 


DEM 


2 


0.053430 


0.0478281 0.048370 1 0.050243 1 0.058673 1 0.065171 I 0.093976 


DEM 


4 


0.057285 


0.0524471 0.053160 1 0 054603 1 0 065050 1 0.072535 1 0.109106 


DEM 


6 


0.060711 


0.057027 1 0.057885 1 0.060060 1 0.070461 1 0.080129 j 0.126889 


DEM 


9 


0.063686 


0.060515 1 0.061369 1 0.064398 1 0.075353 1 0.086016! 0.147006 


DEM 


14 


0.066555 


0.064220 1 0.065201 1 0.068341 j 0.080461 j 0.091432 1 0.161773 


DEM 


19 


0.068085 


0 066183 1 0.067310 1 0.070635 1 0.084379 I 0.096031 1 0.194856 


DEM 


29 


0.071368 


00755781 0.07729910.08172610.104493 1 0.123304 1 0.192463 



Note: forward curves were obtained from the current curves by bootstrapping method. 
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Table of foreign exchange rates (spot and forward): 



Currency 


USD per Ccy, 


USD per 




spot 


Ccy, fwd 1 yr 


ATS 


0.0812843 


0.0811934 


AUD 


0.7000000 


0.6672129 


BEF 


0.0277393 


0.0276916 


CAD 


0,7058159 


0.6984693 


CHF 


0.7041758 


0.7163023 


DEM 


0.5720824 


0.5717753 


DKK 


0.1502765 


0.1495353 


ESP 


0.0067732 


0.0066965 


FIM 


0.1895950 


0.1886958 


FRF 


0.1708701 


0.1706838 


GBP 


0.6898000 


1 .6363808 


HKD 


0.1293494 


0,1221799 


IDR 


0.0002805 


0.0000571 


lEP 


1 .4936000 


1.4107878 


ITL 


0.0005838 


0.0005726 


JPY 


0.0079494 


0.0083655 


MXN 


0.1222046 


0.1055812 


MYR 


0.2928258 


0.2089619 


NLG 


0.5076400 


0.5070257 


NOK 


0.1407856 


0.1371640 


NZD 


0.6266000 


0.5710057 


PTE 


0.0056032 


0.0055445 


SEK 


0.1317020 


0.1275255 


SGD 


0.6327912 


0.5657095 


THB 


0.0252048 


0.0166732 


USD 


1 .0000000 


1 .0000000 


XEU 


1.1342000 


1.1190554 


ZAR 


0.2069536 


0.1866202 
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Covariance matrix of country-industry indices (fragment, actual size 211 by 2 1 1 ): 



Covariance 
matrix 


















Country 


i AU 


1 AU ! AU 1 AU 


Index name 


Country 


Industry 


AU GNRL 


1 AU_BFIN 1 AU BMED 1 AU CSTR 


MSCI Australia 
Index rciAU) 


AU 


AU_GISIRL 


0.0002008 








ASX Banks & 
Finance Index 
(ABII) 


AU 


AU_BFIISI 


0.0001658 


0.000370 

y 






ASX Media 
Index ( AMEI) 


AU 


AU_BMED 


0.0002384 


0.000144 

4 


0.0005618 




ASX Building 
Materials Index 
(ABMI) 


AU 


AU_CSTR 


0.000171 


8.437E05 


0.000192 


0,0004216 


ASX Chemicals 
Index (.ACM) 


AU 


AU^CHEM 


0.0002567 


0.000162 
4 


0.0003524 


0.00023 


ASX Energy 
Index ( AXEY) 


AU 


AU_EIMRG 


0.0001518 


0.000107 

9 


0.0001871 


0.0001241 


ASX Food & 

Household 

Goods 

Index (.AFHI) 


AU 


AU_FOOD 


0.0002046 


0.000192 

3 


0.0002408 


0.0001694 


ASX Insurance 
Index ( All!) 


AU 


AUJNSU 


0.0002008 


0.000195 


0.000255 


0.0001246 


ASX Paper & 
Packaging 
Index 
(.APPI) 


AU 


AU^PAPR 


0.0001811 


7.415E-05 


0.0002756 


0.0002049 


ASX Transport 
Index ( ATII) 


AU 


AU_TRAIVJ 


0.0001754 


0.000138 

9 


0.00025 


0.0001775 


MSCI Austria 
Index (.CIAT) 


AT 


AT_GNRL 


5.911E-05 


9.668E-06 


0.0001002 


7.506E-05 


MSCI Belgium 
Index (.CIBE) 


BE 


BE_GNRL 


8.795E-05 


5.433E-05 


0.000145 


9.354E-05 


MSCI Canada 
Index (.CICA) 


OA 


CA_GNRL 


0.0001152 


0.000125 


0.000127 


7.762E-05 


Toronto SE 
Automobiles & 
Parts Index 
(TlPA) 


OA 


CA_AUTO 


7.868E-05 


8.753E-05 


8.144E-05 


4.668E-05 


Toronto SE 

Financial 

Services 


CA 


CA^BFIN 


0.0001036 


0.000180 

5 


0.0001739 


0.0001015 


Toronto SE 
Broadcasting 
Index (.TCMB) 


OA 


CA_BMED 


4.429E-05 


1 .239E-05 


4.918E-05 


6.247E-05 


Toronto SE 
Cement & 
Concrete 
Index (TIPC) 


CA 


CA_CSTR 


0.0001121 


7.239E-05 


0.0001713 


8.249E-05 


Toronto SE 
Chemicals 
Index 
(TIP2) 


CA 


CA_CHEM 


0.000119 


0.000142 

2 


0.0001386 


8.82E-05 
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Transition probabilities of credit rating migrations (Moody-8 rating system): 



Rating 


Aaa 


Aa 


A 


Baa 


Ba 


B 


Caa 


Default 


Aaa 


0.9338 


0.0594 


0.0064 


0.0000 


0.0002 


0.0000 


0.0000 


0.0002 


Aa 


0.0161 


0.9055 


0.0746 


0.0026 


0.0009 


0.0001 


0.0000 


0.0002 


A 


0.0007 


0.0228 


09242 


0.0463 


0.0045 


0.0012 


0.0001 


0.0002 


Baa 


0.0005 


0.0026 


0.0551 


0.8848 


0.0476 


0.0071 


0.0008 


0.0015 


Ba 


0.0002 


0.0005 


0.0042 


0.0516 


0.8691 


0.0591 


0.0024 


0.0129 


B 


0.0000 


0.0004 


0.0013 


0.0054 


0.0635 


0.8422 


0.0191 


0 0681 


Caa 


0.0000 


0.0000 


0.0000 


0.0062 


0.0205 


0 0408 


0.6919 


0.2406 


Default 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


0.0000 


1 .0000 



Note: Any transition probability matrix may be used (not necessarily 8-state). 
Example 1 

Sample bond portfolio: 



Exposure Name 


Obligor Name 


Rating 


Currency 


Recovery Rate 


912810EN4-11 


AAAAKOT 


Aaa 


USD 


40 


021654878440001 


AAAAAUD 


Aa 


CAD 


50 


19971027-00539 


AAAAKOT 


Aaa 


USD 


40 


027054867160001 


AAAABPD 


Aa 


CAD 


60 


027052668030001 


AAAABPD 


Aa 


CAD 


38 


021655365450002 


AAAAAUD 


Aa 


CAD 


32 


021655355640001 


AAAAAUD 


Aa 


CAD 


17 
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Exposure Name 


RecovRatc 
StdDcv 


Amount 


Matuntx' 


Rate 


Freq- 
uency 


912810EN4-11 


25 


5000000. 00 


15/11/2022 


7 , 625 


2 


021654878440001 


25 


15151000.00 


11/03/2003 


8 


2 


19971027-00539 


25 


14000000, 00 


15/08/2027 


6.375 


2 


027054867160001 


30 


4905000. 00 


15/08/2020 


11 


2 


027052668030001 


23 


1300000.00 


15/10/2004 


11 


2 


021655365450002 


20 


25000000.00 


02/06/2027 


7,6 


2 


021655355640001 


11 


15420000.00 


01/12/2005 


8.25 


2 



"Rate" is the fixed coupon rate of the bond, "Frequency" is coupon payment frequency per 
year. Obligors: 

AAAAKOT - General Electric 
AAAABPD- Hydro-Quebec 
AAAAAUD - Ontario, Province Of 

Sample output (obtained by running CreditVaR program on the above portfolio): 
Portfolio: Exp7 

Number of simulation runs: 100000 
Number of exposures: 7 
Number of obligors: 3 
Evaluation date: 14/04/98 
Started: 11:14:37 
Finished: 11 : 16:24 
Simulation time: 00:01:47 
Analytic solution time: 00:00:00 
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Portfolio Mean (analytic) = 70144493.80 

Portfolio StDev (analytic) = 523028.46 

Portfolio Mean (simulation) = 70143025.60 

Portfolio StDev (simulation) = 530826.76 

Portfolio Percentile (10%) = -105932.78 

Portfolio Percentile (5s) = -297306.87 

Portfolio Percentile (1^) = -812778.79 

Portfolio Percentile (0.5^,) = -970228.53 

Portfolio Percentile (0.1^) = -4912661.85 



Obligor 


Rating 


Present Value 


Mean 
(analytic) 


Mean 

(simulation) 


AAAAKOT 


Aaa 


19190259.30 


19093143.19 


19093014 .25 


AAAAAUD 


Aa 


45880465.72 


44557820.87 


44556843.65 


AAAABPD 


Aa 


6716077.99 


6493529.74 


6493167.70 



Obligor 


StDev 
(analytic) 


St. Dev. 
(simulation) 


Marginal 
StDev 


Marginal 
Risk 


AAAAKOT 


73207.89 


73760.89 


3636.56 


0.019% 


AAAAAUD 


507005. 01 


515008 . 37 


422137.89 


0 . 920% 


AAAABPD 


71014 .55 


78555.92 


12079.87 


0.180% 



Example 2 (actual program output) 
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Exposure data (fragment of input, actual size is 1863 exposures, 207 obligors) 



ID 


Exposure 
Name 


Obligor 
Short 
Name 


Credit 
Rating 


Currency 


Recovery 
Rate 


Recovery 
StDev 


Amount 






rrequency 


15 


15656 


313122 


Aa 


AUD 


40 


25 


1000000 


3/20/08 


6.4 


2 


16 


15657 


313122 


Aa 


AUD 


40 


25 


1000000 


3/20/13 


6.57 
5 


2 


17 


15658 


313122 


Aa 


AUD 


40 


25 


1000000 


3/20/18 


7 


2 


18 


15659 


313122 


Aa 


AUD 


40 


25 


1000000 


3/20/28 


7.05 


2 


19 


15660 


313123 


A 


AUD 


40 


25 


1000000 


3/20/99 


5.41 


2 


20 


15661 


313123 


A 


AUD 


40 


25 


1000000 


3/20/00 


5.41 


2 


21 


15662 


313123 


A 


AUD 


40 


25 


1000000 


3/20/01 


5.64 


2 


22 


15663 


313123 


A 


AUD 


40 


25 


1000000 


3/20/03 


6 


2 


23 


15664 


313123 


A 


AUD 


40 


25 


1000000 


3/20/05 


6.27 


2 


24 


15665 


313123 


A 


AUD 


40 


25 


1000000 


3/20/08 


6.53 


2 


25 


15666 


313123 


A 


AUD 


40 


25 


1000000 


3/20/13 


6.71 


2 


26 


15667 


313123 


A 


AUD 


40 


25 


1000000 


3/20/18 


7.14 


2 


27 


15668 


313123 


A 


AUD 


40 


25 


1000000 


3/20/28 


7.2 


2 


28 


15669 


313124 


Baa 


AUD 


40 


25 


1000000 


3/20/99 


5,61 


2 


29 


15670 


313124 


Baa 


AUD 


40 


25 


1000000 


3/20/00 


5.6 


2 



Example of the program output - Portfolio Summary (analytic engine v^as not used) 
(see Figure 2) 

Example of the program output - Obligor Summary 
(see Figure 3) 

Example of the program output - Distribution of the portfolio forward value 
(See Figure 4) 

Example of the program output - Marginal Risk vs. Exposure 
(See Figure 5) 
Appendix 2 to Part 1 
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Correlation model 
2.1 Market and specific components and country/industry weights 

Let us assume that company 's standardized equity return r of the firm A can be expressed in 
the form: 

f = fii^'-x, + Psl-s > (1) 

where r„ is standardized return of the market component of r and r, is standardized return 
of the idiosyncratic component, r„ and /• are independent. From this follows that 
coefficients /?^, and satisfy condition 

A7+A-=l (2) 
We also assume that return of the market component r„ is linear combination of 
standardized country/industry index returns r, . which are independent of . 

=/?,'; +... + /?,/•„. (3) 

<^or{r,,r^)=0^ i = \ „. 

Equations (I), (3) can be obtained by regressing A equity returns against index returns in a 
usual way. 



ic return . 



There is another way to represent r in terms of index returns r, and idiosyncrat 
This representation is consistent with (l)-(3) but based on countiy/industry weighting 
coefficients. We show below how this representation can be obtained from (l)-(3). 
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Denote by cr, cr, - standard deviations of the firm A equity return and index / return 

respectively, p,^ correlation between index returns r, and ii,J = \ /i). Then from (3) 

one has 



cr. 



o-.r. +... + ■ 



y5„ 



or 



rv/ = + • • • + ^ jCw'.f.r, + . . . + H'„<r,/-, ) . 



(4) 



where >v, , ^^'n are given by 



^=1 // 



(5) 



1/ 



Note, that 

H', + w, + . . . + = L 

Coefficients are called country/industry weights and can be expressed in 

percents. 



From equation (4) follows that 



var 



0, 0 ^ 



n n 



(6) 



Denote 
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(7) 



Then from (6) we have 



Using (I). (4). (8) one can express standardized equity return in terms of country/industry 
weights as 

Hence, to compute this representation one needs to have follow.ng parameters: 
countiy/industry weights, index standard deviat.ons and correlations, specflc coefficient 
(coefficient A, for the marlcet component can be computed from (2) as f},, = VTT^). In 
the current implementation CredhVaR I derives parameters and A, from • specif.c risk 

percent" v defined as a fraction of total equ.ty return changes expla.ned by frm-spec.fic 
movements 



(10) 



2.2 Correlation between two obligors 



Equations („. (3) as well as equation („ ,„„„ ,o con,p„,e ^.^^ 

returns of any two firms. 



Let r<- are standardized equity returns of two different firms B and C respectively and 
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^ J 



S ' 



where 



IE 



=1 / I 



XZv/'u';V,0-^/> V (II) 



A =A/— - A/^^ ^ ' = ^'..,// (12) 
Then correlation between r" and r'* is 

cor(r"y)=f±ft!'P';p, (13) 



Example. 

Suppose we wish to compute correlation between two obligors B and C in the obligor table ir 
Appendix I (Sample Input and Output Data). Obligors shon names are "28664'* and "29715" 
respectively. From the table we have: 

Obligor B ("28664") 

Specific percent .v* = 0.87 

Country/industry weights = 0.5 , u';;,^,.,^ = 0.2 , U',iJ, = 0.3 . 

Obligor C ("29715") 

Specific percent .v*^ = 0.88 
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Country/industry weights w^^ = 0.2 . w^^.b^s- = 0.8 . 



Let us assume following standard deviations and correlations between market indices: 



Index 


Standard 


Correlations 




Deviation 


ID 
GNRL 


ID AUTO 


ID 
FOOD 


US GNRL 


US 
BFIN 


Indonesia, 
General 


0.03 


I 


0.4 


0.2 


0.01 


0.03 


Indonesia, 
Automobiles 


0.02 


0.4 


1 


0.1 


0.02 


0.04 


Indonesia, 
Food 


0.04 


0.2 


0.1 


1 


0.01 


0,02 


US. General 


0.01 


0.01 


0.02 


0 01 


1 1 0.3 


US, Banking 
and Finance 


0.02 


0.03 


0.04 


0.02 


0.3 


I 



According to (7): 

al = 0 03' X 0.5- + 0.02^ x 0.2^ + 0.04' x 0.3^ + 2 x 0.03 x 0.02 x 0.5 x 0.2 x 0.4 
+ 2 X 0.03 x 0.04 X 0.5 x 0.3 x 0.2 + 2 x 0.02 x 0.04 x 0.2 x 0.3 x 0. 1 
= 0.000471 

= VO.000471 = 0.022 

a^c = 0 or X 0.2' +0.02' xO.S" + 2x0 01 x 0.02x0.2x0.8 xO.3 
= 0.000279 

=V0.000279 =0.017 
From (10), (II), (12): 

= Vl-0.87- =0.49. 
/?/=0.87, 

PiD.GNRL " Pks : = yjAy — — = 0.33 , 

0.022 

or» 0.022 
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f>B __ r>B ^^lO, FOOD^ID. FOOD _ in. 0 .> X 0.04 _ ^ ^_ 

PiD.FooD " Pm : u.**y — 7~rrr — - u.z/ , 

cr^ 0.022 



/?;; = Vl-0.88' = 0.47 , 
P'^ = 0.88 , 



A^kav;. = = 0^7^^^ = 0.06 

ar 0.017 



;?c' _ r.c ^4s.8rix^t:s.BFiy ^ .-OSx 0.02 ^ 

A/5./?/./.v = Pa/ = ^ = 0.44 

cr^ 0.017 



Finally, compute correlation between and using (13): 

cor{r^\/ ) = 0.33 x 0 06 x 0 0 1 + 0 33 x 0 44 x 0 03 4- 0 09 x 0.06 x 0.02 + 0.09 x 0.44 x 0.04 
+ 2 X 0.27 X 0.06 x 0.0 1 + 2 x 0.27 x 0.44 x 0.02 
= 0.017 



2.3 Simulation of correlated equity returns 



Obligors correlation matrix is used in the model to generate correlated samples of normally 
distributed random vectors of equity returns. In this section we describe the procedure based 
on Cholesky decomposition of the correlation matrix of index returns 



Suppose there are K obligors and their standardized equity returns are represented as (see 
(12) above): 

Denote by 



(14) 
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Pz\ PZ2 
^PnX P„l 



Px. 



c - 



PnnJ„^„ 



(C„ 0 







P: . 


- pr 




[P's 


0 . 


.. 0' 


B = 


Pz 


P2 ■ 


.. pr 


, A = 


0 


PI . 


.. 0 






P- ■ 


■ p^. 




.0 


0 . 


•• PsL 



where (\, is correlation matrix between indices and is KxK identity matrix. 



Then correlation matrix C of equity returns equals to. 

r = J:{rr')^ liw'rr'w)=W'E(rr')^ = W'C\,IV (15) 
Let us assume that matrix //„ is obtained by Cholesky decomposition of the matrix 

Q/ = ^.\t^lt (16) 

Denote by 



^ 0 / 



(17) 



(// is /:x(// + A') matrix). Then 



C = W'H„Hl,W = HH"^ 
are Cholesky decompositions of Q, and C respectively. 



(18) 



SUBSTITUTE SHEET (RULE 26) 



wo 01/52121 



PCT/CAOO/01440 



38 



For simulation of a sample of K dimensional normal vectors R with zero vector of means 
and covariance matrix C (standard deviations of coordinates of R equal 1) we first simulate 
/J + AT dimensional vector of standard independent normal variables 



From (18) follows that R has distribution iV(0,C): 

cov{R)=j:{Hxx'fr)= HE(xx^)h' = HI „^^H^ = HI^^^H"^ = C . 
Equation (19) can be used for simulating correlated sample of equity returns in the model. 

Appendix 3 to Part 1 

Pricing Procedures 
3.1 Pricing Procedure for Floating Rate Notes 

Floating rate notes are designed so that the coupon payments depend on some current interest 
rate index (usually LIBOR or some other known rate), which is called reference rate. This 
causes the value of a Floating Rate Note (FRN) to be close to par at all times Credit risky 
FRNs pay coupons combined of reference rate plus seme spread over it, which depends on 
the credit rating of the issuer. 




and then compute R as 



R=Hx = W^H.jX 

=(b' a: 




SUBSTITUTE SHEET (RULE 26) 



wo 01/52121 

The notation is as follows: 



39 



PCT/CAOO/01440 



P - price of note, 

r, - expected future reference rate 

V, - spot long-term yields, used to discount the cash-flows at time 

v - stated credit spread, fixed under terms of the FRN 
m - maturity of note in years 
// - coupon payment frequency 

Assuming flat structure of credit spreads, the price of the note per 100 of face value, is 

Here, r, is forward reference rate for the period between /-th and 0 ■ I)-si coupon payments. 
Assuming the reference rate is annually compounded, the way to compute forward reference 
rate n is the following: 

r - ■ - 1 

where R, is the rale on the reference cuive for time point /, . 

Within the current framework of CreditVaR I model, we will take existing current and 
forward corporate curves as a proxy forj/,, since they incorporate market required spread. The 
reference rate R and credit spread .v should be supplied as parameters of the FRN. 
Note that for computing present value of an FRN we use the supplied reference curve, but to 
compute forward I -year value of an FRN, a forward reference curve has to be obtained. The 
rates n in the computation of forward value of an FRN are therefore "forward forward'* rates. 
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3.2 Pricing Procedure for Interest Rate Swaps 

Swaps can be characterized as the difference between two bonds, which can be either fixed or 
floating rate. In the swap agreement, principal is not exchanged, which makes swap different 
from a bond from the point of view of credit risk. Another difference is that in the event of 
default by the counterparty, the sign of the value of a swap determines weather the other 
party pays or receives some recovery amount (this depends on netting rules). 

Assumptions: 

1) There is no more than one credit event per year 

2) The credit event happens one year from evaluation date 

3) Netting is done for positions with the same obligor and seniority 

4) One recovery rate is applied to a netted position 

5) Implementation allows different recovery models depending on the sign of exposure to a 
bank, including recovery rate 0 if the value of the swap to the bank is negative 

6) A bank's credit events are not modelled (i.e., the bank's credit rating is constant). 

The notation is as follows; 

K- value of the swap to the bank, 

5/ - value of the bond modelling the cash flow received by the bank, 

Bj - value of the bond modelling the cash flow paid by the bank, 

r/ - recovery rate applied if bank receives in the event of default by counterparty, 

rj - recovery rate applied if bank pays in the event of default by counterparty, 
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Q - notional pincipal in swap agreement. 
The value of the swap to the bank then is; 

B, and B2 are present (or forward 1 year) values of the underlying bonds, obtained by 
standard procedure for fixed or floating rate bonds. Ideally, swap bears credit risk only if the 
value of the swap to the financial institution is positive. In practice, however, if counterparty 
defaults, it might not free the bank from liability to that counterparty (due to netting rules). 
Therefore, in case of default, the recovery amount is calculated as follows: 

\fV ■ 0, the recovered amount is r,*V, 
If V 0, the amount paid is r;*y. 

Here, we don't take into consideration the event of default or credit rating migration by the 
bank. Recovery rates r, and r, are determined from legal considerations at the time of default. 
We can use Beta distribution again to simulate recovery rates r, and (in current 
implementation, we take n-r,). Alternatively, r, can be put to 0 if the swaps with negative 
exposure are excluded from the computation. 

3.3 Pricing Procedure for Loan Commitments 



A loan commitment is composed of a used portion and an unused portion. Used portion has a 
Loan (drawn fiinds) and a Letter of Credit (undrawn fijncis) parts. Interest (usually LIBOR 
plus spread) is paid on the Loan portion of a loan commitm :nt, LIBOR plus LC fee is paid on 
the Letter of Credit part (throughout its lifetime), and an Urused Fee is paid on the unused 
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portion of the loan connmitmem. There is also a Facility Fee applied to the total loan 
commitment amount throughout its lifetime. A 1-time Upfront Fee paid on the total loan 
commitment amount at the initiation can be ignored in the Credit VaR I framework since it 
bears no credit risk. 

The amount drawn at the risk horizon is closely related to the credit rating of the obligor If 
the obligor deteriorates, it is more likely to withdraw additional funds. On the other hand, if 
its prospects improve, it is unlikely to need the extra borrowings. 



Let /'be the total (authorized) loan commitment amount: /- - the amount drawn (in loans); LC 
- the amount in Letters of Credit; II - unused portion of the loan commitment. 
We have thai 

Let / be LIBOR and .v be spread over LIBOR, so (I - s) is the interest paid on the loan; Ic - the 
fee paid on the letter of credit, uf- the fee paid on the unused portion;j^- facility fee applied 
to the total authorized amount of loan commitment. Then present value (in the CreditVaR 
framework) is computed as follows: we treat the loan and letter of credit parts as floating rate 
notes (because of the way interest is paid); facility fee and unused portion fee are valued as 
bonds with no principal repayment. 



PV ^ 



'^[/ X uf / n 



^ LC X Ic / n 



+ 



Txf//n 



where y is the annualized rate taken from the corresponding current corporate curve used to 
discount cash flows (depends on credit rating of the obligor in the CreditVaR framework); n 
is the frequency of coupon payments, and m is the number of years to maturity. Forward 
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value (FV ) is computed similarly, with y being the forward (1 year from now) rates taken 
from the corresponding corporate curve. 



If downgrade or default happen to a particular obligor, it is likely to withdraw additional 
funds from the loan commitment line. To model this we assume that there is an "expected" 
drawdown (AT)), or "average commitment usage" if the obligor changes credit rating (see 
Credit Metrics Technical Document, p.45 for appropriate information): 



Credit rating 


Averace commitment usace (ED) 


Aaa 


0.1% 


Aa 


1 6% 


A 


4.6% 


Baa 


20.0% 


Ba 


46 8% 


B 


63.7% 


Caa 


75.0% 



Then in the case of any credit event (downgrade, upgrade, default) the used portion of a loan 
commitment is adjusted in the following way: 



Downgrade 

Let XVo (y% ) is average commitment usage given in the Table above for the current (new) 
credit rating respectively, X < Y . Then new unused, loan and letter of credit portions are 
computed as 

\-X ' 

- new used amount T-U^ is proportionally distributed between loan amount L^^ and 
letters of credit amount LC 
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Upgrade 

In this case X >y . New unused, loan and letter of credit portions are computed as 

A_ = max{/.-(f/„^-f/)^0) 
/'C_ = /'-/-_.-^/„... 

Default 

In the case of default it is assumed that used ponion is the maximum of current used portion 
and 90% of the authorized amount 7" : 
= max(/. + /,C,0.9 x '/ ) . 

Forward value of the loan commitment in the case of upgrade or downgrade is computed as 

In the case of default the loan commitment value is x R , where R is the recovery rate 

sampled from the Beta distribution (consistent with CreditVaR I framework), and total 
forward value is 



SUBSTITUTE SHEET (RULE 26) 



wo 01/52121 PCT/CAOO/01440 

45 

Part 2 - CreditVaR Methodology: Credit Migration Process Under Two Probability 
Measures 



As described in Part 1. the CreditVaR is a model for measuring and analyzing credit risk in a 
portfolio context. The methodology used is based on credit migration analysis and on the 
assumption that an obligor's credit migration is driven by its asset value. Correlations 
between equity returns are used to compute joint probability distribution of obligors' credit 
migrations. Each obligor s standardized equity returns are decomposed into weighted average 
of market indices returns (multi-beta model) where the weights are specified to appropriately 
reflect the obligor s participation in the corresponding markets and to model obligor's 
idiosyncratic returns. 

In the Monte Carlo simulation version the model calculates distribution of the portfolio 
values, statistical characteristics and risk measures defined as percentiles at given confidence 
levels (95%, 99%. etc). Simulation engine generates scenarios based on the "real world" 
probability distribution (P) of the risk factors. The distribution is assumed to be lognormal 
with parameters computed from the historical time series data. For each generated scenario 
the ponfolio value is computed under the martingale probability measure (Q) The 
relationship between these two probability measures is derived from the assumptions about 
stochastic processes for index and stock returns. In particular, we derive the relationship 
between P and Q probabilities for obligors to migrate from one credit class to another 
including probabilities of default (transition matrices). This allows one to analyze portfolios 
of different instruments including credit derivatives. 
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In the Credit VaR framework it is assumed that credit migration process forms a discrete 
Markov chain with fixed moments of time T = < <..,<(/, <...}, fJ^ ^ 0 ■ In sections 
1-5 of this Part we present mathematical models of the migration process and define 
probabilities under P and (J measures to migrate from one credit class at time to another 
credit class at time Theoretical results of sections 1-5 are applied in sections 6-8 to 

compute P and O transition probabilities from market credit spreads and credit migration 
historical data. 

1. Stochastic processes for market index and stock returns under F measure 

Suppose that IS -{B\li' B^^) is a standard /^-dimensional Brownian motion on a 

probability space (n,F,/^) {B\B^' B^^ are independent), where F = {/-^,/>0} is a 

standard filtration of B , 

Let I = {/',/- y } be a set of all market indices. Denote by the index I'" value at 

time / > 0. It is assumed that I]"* satisfies stochastic differential equation (SDE) of the form: 

(1.1) 

^7 = v-, m = \0 M . 

where //f , 0^^ are piecewise continuous deterministic functions of / , / = 1,2,...,ZJ) . 
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Let us fix some time interval ',] corresponding to two consecutive time steps from f : 
= f/; > = f/,., , / > 0 . Denote = Ir(F,")- ln(F,;) be the index /" log-return for the time 
interval [/„,/], t„<t<.t^. Then from Ito's lemma and (1.1) follows that 



(1.2) 



for /„ < / < /, with initial condition jt,"^ = 0 . In a matrix form (1.2) can be written as 

t 

= + \&MR, ./.,:£/</,. jc, = 0 . (13) 

where 

^.=(*; ^f. 





dB, = (ub: JB^J 

Let us consider obligors (issuers) S = tv^^^...,^-}. and let //,- be an obligor S- stock 
price at time / We assume that //," satisfies SDE: 

dH: = H:[a:dt + ^.dl?; + al^dS^ +... + al^dB^), , > o. 

(1.4) 
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(1.5) 



where a" , a^*, are piecewise continuous deterministic functions of time, / = 1,2 D , Then 

for the stock log-return = ln(//; )- ln{//;^ ) one obtains 

oL - 

0 

^o^/^/,. =0. 

Similar to (1.3) this can be written as: 



(1.6) 



where 

y, = {y) >f 



2. Multi-beta model 

In the CreditVaR framework correlations between stock returns are derived from the index 
return correlations using multi-beta model. This is linear regression model that represents 
stock returns as linear combinations of index returns and stock specific (idiosyncratic) 
components: 



(2.1) 
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« = 1,2 A'. 

where e" and .v" are uncorreiated: co\{s",x"')=0 , and E{e")=0, m = \,2 M, 

n = 1.2 . In the matrix form: 

(2.2) 

where A" = (A'" A^^f, A =||A"-1L.,, is a matrix of stock "betas", =(^; s^'J is 

a vector of idiosyncratic components which are uncorreiated with index returns: 

cov(£,,.v,) = ||cov(<,.v;'|^.^^ =0 

(2.3) 

(here and below expected value E(.) and covariance cov(.) are conditional to cr -algebra 

It is assumed that rank(d^)-M <D for all (Note, that if rank{o,)<M <D for 

some / = / then it can be shown that there exists time interval [^„./,], / e [f.?] and index 

subset I el, 1^0, such that for all indices in I their log-returns can be represented as 
linear fijnctions of index I.I log-returns with time dependent deterministic coefficients. 
Then index set I can be reduced to subset III on [/".?] ) Then from (1.8), (19) it follows 
that the (A^/ xA^) -matrix OfiJ is nonsingular on and 



V'o A'« ) 



(2.4) 

Indeed, from (1.3,(1.6) follows that 
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cov{x„x,)= jO^eJc/s, 



(2.5) 



(2.6) 

Then (2.2) and (2. 5) give 

coxiy; .x,) = cov(pl' + p,x, + £,.x,)= fi, cr»<.r, . .r, ) = fi, . 

'» 

Therefore, 

] a, Olds = p,\efilJs 

'» In 

and (2.4) holds. 

Equations (2.2) and (2.4) also give the covariance matrix between equity specific 
components: 

cov{£„£,) = cov{y, - /?;• - /}c,.y, - P: - fix, ) 



(' Y' VY' 

= \o,c7\ds - \afilds \efilds jO,a[ds 



(2.7) 



In the special case of constant volatility matrices 0 and a- one has 

where is (d x d) identity matrix. In this case equity specific components are 
uncorrelated if and only if a[/^ -o'^{o 0^ X o)j = , where X is a diagonal matrix. This 
means that orthogonal projections of the rows of stock volatility matrix a on the linear 
subspace, which is an orthogonal compliment in of the linear space generated by rows of 
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O, are orthogonal. Note, that in th s case covariance matrix of stock returns equals 
<=o<\ ' y, ) = (' - 'o iAoeffi"^ + a). 

3. Credit migration model under probability measure P 

It is assumed that at any time ^/ e f each obligor S" e S belongs to one of credit classes 

^^=^'^1 ^f:)^ where corresponds to default. Suppose that credit migration process 

forms a Markov chain on R with discrete times f and absorbing state . Denote by R,". a 
credit class of obligor .V" at time U ef . Similar to section 1, we fix a lime interval [/„,/,] 
corresponding to two consecutive time steps from f : /., =(/, . /, =f/,., . / > 0. Suppose that 
K = P'j {'>• -> ) is a probability of S" to migrate to credit class R^ at time /, . 

(3.1) 

Denote by 

(3.2) 

a transition matrix of S" for the time interval [/„,/,]. Then 

Xp"M>^0=i y = i K 

• /'L = 0. 7 = 1 K-]. 

For each nondefault state /< and probabilities p^Jf,, ->/,), y = l K define thresholds 

Kiif»A)t j = l,...,Ar - 1 , for the standard normal distribution as follows: 
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(3.3) 

where N{.) is standard normal cumulative density function. 



Let »»■, = Z ;.',{}', - /■-'(>',)) be a vector of standardized stock returns with correlation matrix 
given by 

f, •) 

cor(»;.M;) = Z;,cov(y,.yjZ:] = Z;', Jo-^crV.v Z.. 



where Z,. , is a {N x N) diagonal matrix with stock returns standard deviations 

7var(v,' ) 0 ••• 0 



0 Vvartv;) -- 0 
0 0 — Vvarlv',^)J 



(var(.) is a variance under measure conditional to o" -algebra ). 



(3.4) 



(3.5) 



Similar to the CreditMetrics framework, the Credit VaR credit migration model is based on 
the assumption that credit class of the obligor S" at time /, is defined by the value of its 
normalized stock return u," and thresholds //."^(/o-'i) 



SUBSTITUTE SHEET (RULE 26) 



wo 01/52121 PCT/CAOO/01440 

Assumpiion: Given R^^^R^, R;; = if and only if e fc, X A, V,(/o,Oj, where 
'J^h2 = 00, A;^, (/,,/,) = -00. 

From the Assumption it follows that joint migration process for obligors can be defined using 
a joint probability distribution of their normalized stock returns. Indeed, suppose that 

K = ^^f, ' e S , // = I . Then normalized obligors' stock returns have 

dimensional normal distribution with zero means, standard deviations equal to 1 and 

correlation matrix (4 4) Denote by /^(x,,x. x,.) their joint density function. Then the 

joint probability that at time /, obligors will be in credit classes respectively 
is given by 

^(K-^. <=/?,J/i;>/<. < = «,J 

(3.6) 

4. Stochastic processes for index and stock returns under Q measure 



In this section we derive stochastic differential equations for the index and stock 
under the martingale probability measure O. 



returns 



Let r, be an instantaneous risk-free interest rate at time /. We assume that is a 

deterministic positive function of /. Then index values f;", /w= 1,2 A/, and stock prices 

//", // = 1,2,...,// , for non-dividend-paying stocks (this assumption can be relaxed) satisfy 
following equations: 
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dvr = 17 (/■ Ji + eyh] + 0^\JB- +...+ e:;^JB^ ) , / > 0 . 

(4.1) 

I ;."■ = »•"', /w = i.2 .V/, 

JH: = H:[r,ih + c7;,dH]+a\JH; +...-i-cj'^^^^^ />0, 

(4.2) 

= /, = 1.2 A', 



where 

= /y/-rJr;'c/v, d=\....j). (4,3) 

n 

are Brownian motions under the probability measure O defined by the Radon-Nikodym 
derivative. 



^ = expf-Xjr;'^/^:'-^l|r,r./vl 

V ''■111 -n J 



(4.4) 



where vector-function t,={t],t; r,"f solves equations: 

C^-; +^0 r + ... + B-^zl' = //,"' - /; , m = 1,2 M , 



+ - . . . + <„r," = a," -/;,//= 1,2 N . 



(4.5) 



(4.6) 



In the matrix form: 
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(4.7) 



where fi,={^\-r, M^'-rJ, «,=(«,'-'•,. 

(4 7) exists for all / e {/„./,]. 



-r,J . We assume that solution r, of 



Similar to (1.5), (1.8) one computes index and stock log-returns from (4 1), (4.2) for the time 
interval [foJ^ = [U,M,,,]. 

t 



(4.8) 



(4.9) 



where 



=[\[r.-{ioijAo\j-.Mo\j)]is ^[^4k;^k'J.....(^-y)}/.vJ. 

*.=(*; A.^r 

<JB,={fiBl dBf] . 



Let E^{) cov'^{.,) be an expected value and covariance matrix under the probability measure 
O, conditional to cr -algebra F,^ . Then from (4.8), (4.9) follows that 
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E^'{x,) = a, (4.10) 

Fy{y.) = h, (4.11) 

cov^'{x,.x,) = cov{x,,x,)^ \d^el Js, (4. 12) 

f 

cov^{}\.x,)^cov{y^.x,)= ja^d^ds , (4.13) 

f 

cov^'{y,.y,) = cov(jr,y)^jcr^cylds'. (4.14) 

In particular, 

var^^(>;)= var(>;) (4.15) 

for all // = 1,2 jW . Moreover, 

COV^ (f, . X J = CYM<^, . -V, ) = 0 . (4.16) 

c«v^^(6; . 6; ) = cov(^, . f J . (4.17) 

but 

K^'{i:,)=h,-/3/i,-/3';^0. (4.18) 



5. Credit migration model under Q measure 

In section 3 we defined the relationship between credit migration events and movements of 
the normalized stock returns. This allows to compute probabilities of credit migrations for 
different credit classes under the (J measure. 
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Let obligor e S be in the credit class R;^ = /< at time . According to the assumption in 
section 3 its credit rating at time /, will be R^^ = R^ if and only if >v,: e fc(r„,/,).AV,(/,./,)J. 

^ = thresholds /»:,(/„,/,) are derived from the standard normal distribution. 

Then probability under O measure for obligors S\S'- to migrate from the credit 

classes /<_ = /<_ Rj^ = /^^ at time /, to the credit classes 7^; = R^^ R^^ = r^^ at time /, 

respectively is 

<=/^,j/<=/<. 

= <M ^ k., (A../. ).//.;.„ , (a,,/. ))..., e fe.^^. (,,..,,)]) 

(5.1) 

As follows from the definition of = (»■,' w^'J and (4 15) 

= »';+c:(/„,/) 



(5.2) 



where 



(5.3) 



Q('o.O is a constant given by 
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[» 

I 

Then for the vector = n-^ ] using (3.4), (4. 14), (4. 15) one has 

COV^^(u*, , ir, ) = cov^^ (u; , vt; ) 

= z;'cov(>',>',)z;.', 

= COv(u',.»f, ). 

(5.4) 

From equations (4.9), (5.3). (5.4) it follows that the random vector h-, =(»I-' w^J has a 

normal distribution under the probability measure (J with zero vector of means, covariance 

matrix cov'^l*!', ,u-, )= cov(»i', ,m-, ) and probability density function /^.(.v,,.r x^.). Then 

w" is a standard normal variable for each // = I,...,.V and 

<M. ^o=cX»-;: ete(A..o.Av,(A,./.)j) 

=c?k e {KM-c:,{tMKj ,('o.O-c(a../.)]) 
= Ma".-, ('...', )- ))- A^fe (/„./, ) - c,:(/.,./, )). 

Therefore, similar to transition probabilities p"j{t„ — */,), transition probabilities (/"j{f„ — 
under martingale measure O are defmed by shifted thresholds 
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8:.M,.|^ ) = K. ) r 

(5.6) 

Note that transition probabiiities -> 'i) depend on // and can be different for different 

obligors. Equations (5.5), (5.6) show the relationship that exists between two transition 
matrices 'C ('., ) = Ip.:, (a. and = |./:, (/., -^/,|. 



Similar to (3.6) for the O probability of joint migrations we have 



6. Construction of transition matrix under 6> -measure over risk horizon 



In this section we suggest, as an application of the theory developed in earlier sections, a 
practical way to obtain a Q transition matrix for the period [/„./,], t^=i^^ + H , where H 
denotes a risk horizon, typically 1 year. We suppose that we are provided with a P transition 
matrix and a set of Q default probabilities, one for each initial credit rating, over the period. 
The matrix obtained this way will be used in subsequent sections. 
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The main point here is that of implementation issue. In practice, the differences between 
default thresholds of the given P transition matrix and those determined by the given O 
default probabilities may not be independent of the credit ratings, contrary to the theory (see 
(5.6)). This is the case usually because the P transition matrix is obtained empirically from 
past historical data while the O default probabilities are obtained from current price data. 
Below, we make a modification to the given P transition matrix and then apply the theory to 
obtain a (J transition matrix over the risk horizon Note that the last column of such O 
transition matrix will coincide with the given O default probabilities. 

First, we introduce some general notation. For obligor .S" € S and observation times / >.v, 
.V,/ G T , let q"{s ^1) denote the (J default probability of .V" over the period [sj] given its 
initial rating of at time s , i.e., 

v:cv ->/)== q{r: = D I r: = /<) , / = 1,2 k . 

Recall that the default state 7^^. = D is absorbing so that ^/^.(.v /) = 1 . Put k = AT - 1 , and 
package them in a k -dimensional vector 

we have q''Cv / ) > 0, where 0 denotes the k -dimensional zero vector. 
Now suppose we are given a P transition matrix T^P =^ \- p^^ \\ for each individual obligor 
S\ respectively, and (J default probabilities q''(/„ ->/,) over the risk interval [/,„/,]. Put 
= 7: /, ) and r; = (^f ) , 1 < / < k . Let 
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c: = ( -.y/.. ; = 



[l((<J-'-k:y-^...-(a:.,r)/.vj-' 

as in Section 5. Define transition matrix /„'' = || p"^ {i„j, ) || by 

where u, = '' J " ; y,;;;.^ (/,.,,, ) = o if 7 < A' and /;;;,,. (/..,/, ) = 1 . We have adjusted 

the last column of the given matrix T/: so that the new P default thresholds satisfy 
'(/-^/"/.(^..^ 7 for / < A\ i.e., they differ from the (J default thresholds z" by 

constant C." . The other entries corresponding to the probabilities of migration into non- 
default states are then adjusted while maintaining, for each row, the same relative weights 
with respect to the survival probability, i.e., 

■-/^V(/„.',)=l^ '-/<^ (6.2) 

Now, we apply the theory of change of measure to /;,'' to obtain the desired (J transition 
matrix over the risk horizon. 

With 7;; = (I p:^(t,„t,)\\ of (6 1). let //.".(A,./,) be its thresholds. In particular, we have 

KAf..^',) = A^"'(/<\(A„/,)) = --" + (.'.:. (6.3) 
Let 7,y = II c/"^ ) II be the corresponding transition matrix under (J ; by (5.5) and (5.6), its 
(/,./) -th entries are defined by 

('.„',) = A^a-:, , c. ,/,))- a^(^:, (/„ ,/,)), /. =1,2 k (6.4) 
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where ^f,"^ (/.,./, ) = /i,"^ (/„./, ) ) . Using (6.3), (6.4) and N '(—x:) = 0 , we have 

<i:.K ('.. ^f^) = f^(K< (/.„/,)- f ',r ) - ^(-x ) = N(z: ) = (/., ^ > 

for each / = 1,2 k . Thus, by construction, the last column of /'^ fits exactly the given 

default probabilities under (J over the period [/„,/, ] : 



(/.../. 


















'■/;./.(/„.^) 








[ 1 J 


I ' J 




1 





(6.5) 



In summary, a suitable adjustment is made to the given /* transition matrix such that its last 
column is compatible with the given (J default probabilities under the theory, i.e., their 
thresholds differ by a constant T independent of initial credit rating class. For each row, the 
remaining entries in the transition matrix are also adjusted by retaining same relative 
weights with respect to the survival probability. The constant (' is then subtracted from other 
non-default thresholds of the modified matrix to obtain the corresponding thresholds, and 
therefore the entries, of the desired (J matrix. 

Note, that in most practical applications transition matrices for individual obligors are not 
available. Instead only one common for all obligors P transition matrix is given, it can be, for 
instance, transition matrix published by one of rating agencies like Moody's or S&P, or the 
internal bank's transition matrix. For such case following modification of the procedure 
above can be suggested. We assume that (J default probabilities q"(/„ /, ) are also obligor 
independent so that q"(^„ = q(/., ->/.). Let be the / -th component of q(/., /, ) , and 
put 
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<^'o=-|:i:(n^,-'). 



Here, for ('„, we may want to average the differences w, -z, over only the / -th terms for 
which 7, and p, ^ are not too small (so that the thresholds are away from - oc); also, instead 
of taking the average, one can choose C„ which gives the least mean square error. Similar to 
(6. 1 ), define the new transition matrix T'' = || ^ (/„,/, ) || by 

, l-N(z +C„ ) 
where u. = —p . /v,//.,./,) = 0 if j < K and /^,,- (/„,/,) = 1 . Therefore, the last 

column of the given matrix 7). is adjusted in such way that the new P default thresholds 
satisfy A'"' (a.^ (/„./,)) = r, + C„ . i.e., they differ from the O default thresholds r, by 
constant C\ independent of / . Similar to (6.2) we have 

-. ; - = - — r — for all i,j<K. 

^-P,j:i'n>t,) ^-P.K 
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7. Forward transition matrices under O 

In this section, we suggest an algorithm that builds a series of transition matrices matching 
the given term structure of default probabilities over various maturities. Throughout this 
section, we assume that probabilities and transition matrices are under O -measure. Such 
series will be useful, for example, when pricing credit derivatives via risk neutral valuation. 

Recall that for each obligor .V" € S we assume that its credit migration process observed at 

times in T forms a Markov chain on states R under 6) -measure. We view T as 'universal' 

in that all observation times considered from now on are elements of T . In practice, we have 
a particular finite subset of observation times known from beginning. 

Let T = { ;/„ < //, <iia \ (c T) be a set of observation times. By a series of transition 
matrices on T , we mean a sequence of transition matrices over the periods 

K,'/,], ■ , 

With such series, the transition matrix over any period [//^,»//^] with a<h is given by the 
product of transition matrices over ], [/^.,,/^]. respecting the order. 

^et T„ = {/, =/„ +///]/ = 0,1,2 L } be a set of 'standard' observation times, where H is the 

risk horizon; we assume that T„ q T . For obligor .V" g S , suppose we are given a term 
structure of default probabilities under O -measure. By appropriately interpolating if 
necessary, this means that we are given the O default probabilities satisfying 

0<q''(/.,->/,)<q''(/„ -^/,)<...<q'»(/„ -^/,)<e, (7.1) 
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where e is the k -dimensional vector whose components are all equal to i. For /> , 0 < 6 S Z. , 
and an invertible transition matrix A4 over the interval [/..^/^J define vectors ql,(/f, -^ft) 
follows: 



I 



(7.2a) 



for /e {A,A + l,...,/.}. Potentially, these vectors may have negative components because 
A/ ' can have negative entries. However, if these vectors satisfy 

0 = q.u(^ -^/J<q;,(^ )<-.-<q;,(/, -»/J<e, (7.2b) 

then we refer to them as the /ime (Jonx ard) default prahahiliues obtained fnmi the initial 
term structure of default probabilities (7 1 ) usit}^ transition matrix A/ . Note that the equality 
in the first relation of (7 2b) is equivalent to the condition that the first k components of the 
last column of M equals q'*(/„ — >/ft), i.e., it matches the initially given default probabilities 
for the period [/„,//,]. 



We wish to construct certain series of (J transition matrices 

'C'(/n -^/,), n\t, , ^/J (7.3) 

on T„ such that the following proposition is satisfied with Zj - T^\t^^^ —> ^i) 



Proposition: Given initially a term structure of time /„ default probabilities satisfying (7.1), 
there is a set of invertible matrices Z, , Z, , . Z; such that for each A , 0 < A < A , the time 
(f, forward default probabilities obtained from (7.1) using M = Z^Z, •••Z^ satisfy (7.2b), i.e., 
they are non-negative and increasing. 
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Proof of Proposition, We shall produce Z/s recursively. Ai each stage, we aim to keep the 
entries corresponding to the probabilities of non-default migrations as 'close', in terms of 
maintaining their relative weights (This can be replaced by any other reasonable preference. 
For example, instead of seeking to retain same relative weights amongst the non-default 
probabilities, one can skew them in the direction of default probabilities: if default probability 
has increased then one can increase the relative weights towards lower credit classes, etc.), to 
those of the transition matrix which was produced from the given data for [/„,/,] 
Section 6. 



in 



For the initial interval [/.,,/,], let Y\ = ; its last column. 



by (6.7), matches 



the given default probabilities for the interval. If Y, is invenible and condition (7.2b) already 
holds for A/ = , then we set T^\f^, = More generally, consider the following 1- 
parameter family of matrices 

Z,{a) = a\\+(\-a))\, 0<a<l, (7.4a) 

where 



diag(/-;, •-./;') q"(/„->/,) 
0 I 



0 I 



(7.4b) 



with CI, =q:(t„ ^/,)<1 and /-/ = !-<(/., -)■/,) for / = 1,2 k . The matrices Z,(a) ail 

have the same last column and in panicular is independent of a and equal to that of Y\ 



Choose the "smallest" a . call it a,, in (7 4a) such that vectcrs q;,(/, -^/,) defmed by (7.2a) 
with invertible A/ = Z,(a) satisfy (7.2b) with h = \ there (in practice, we fix a positive 
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integer / > i from beginning, and r ;strict values of a in the set a(/) = {0.1.7.2. /.....^/- o./j}. 
Then «, is chosen to be the smallest a in a(/) for which condition (7.2b) holds.) That such 
or, exists is clear since for a = 1 , the condition is satisfied. Thus, we have the time forward 
default probabilities 

0£q"(/, ->/,)^q"(/, ^/,)<...<q"(/, /, ) < e (7.5) 
with respect to M = Z, (a, ) ; here, subscript A/ is omitted. We set 7j'(/„ ->/,) = Z, (a, ) . 

Now, for subsequent periods [/„/,.,]. / > 1 , we proceed recursively to construct invertible 
transition matrices. Using the time /, forward default probabilities (7.5) as given, we can 
obtain in a similar way an invertible transition matrix Z, on [/,,/,] whose last column is 
compatible with the (forward) default probability q"(/, ^ /,) and whose other emries have, 
in each row, the same relative weights as those of T,^' (or, as before, have the weights 
distributed according to some preference) Note that the time /, forward default probabilities 
obtained from the time /, default probabilities (7.5) using Z, of [/,,/,] will be the same as 
the ones obtained from the time A, default probabilities (7 1) using Z,Z, on [/..,/,]. 
Repeating the procedure for intervals [/,,/,], [/,./J and so on. we construct a desired series 
of transition matrices, (7.3). on T„ . 

The key point about the construction is that the series of transition matrices built up this way 
fits the given term structure of default probabilities over each maturity /, , 0 < / < Z. . in that 
the cumulative transition matrix Z,Z, -Z, on the interval [/,„/,] has (the first fc 
components oO the last column matching the given q''(/„ ->/,). 
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Remark. There are many choices for such series. We have suggested simply one method that 
is practical and reasonable in incorporating much information about the P transition matrix 
given for the initial period as well as preserving non-negativity and increasing of the forward 
default probabilities defined by the constructed transition matrices. Other alternate methods 
for producing such series can be made available. Actually, in practice, one may prefer to be 
given as exogenous input a series of forward O transition matrices in place of the initial term 
structure of (J default probabilities. 



8. Joint process under O on T and application to pricing of credit swaps 

In the previous section, we have built series of O transition matrices on standard observation 
times T„ for single processes R" . In application, for example in pricing of credit derivatives, 
we need to 'translate' this series to cover a different set of observation limes 
T = {//,,<«,<•••<: w„} (T includes nonstandard' time steps corresponding to payment 
dates, times to instrument s maturity, etc., which need not belong to T„) for a joint process. 
For this, it is sufTicient to describe it for single processes, for then the joint migration Markov 
process will be determined by (5.7) on each period ] as the co-migraiion of the 

obligors are governed by the correlation among their asset returns. We make the assumption, 
in lack of any better information, that this correlation is constant in all periods, but this too 
can be given as input and made different for different periods. For the single process, it 
essentially requires that we provide a way to 'interpolate' a transition matrix for a given 
period to define a transition matrix for a subinterval of that period. We shall simply take a 
certain Taylor series expansion of the 'fractional' power of the transition matrix, the power 
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being proportional to the subintervars length, and when necessary make a minor adjustment 
in some cases. 



Let Z be a transition matrix for the period [/^,/.,., ] and suppose [//, v] c [/^,/.,., ] is a proper 
subinterval; put 9 = Now consider the matrix Z" defined by 

=i;((e)),(Z-/)* , (8.1) 

k'tt 

where ((6))^ =0(0- l)---(0-/: -hi)/*!, as in the Taylor expansion around the identity matrix 
/ When scries (8 1) converges, we may truncate it after a certain number of terms to 
evaluate /" with a small error. Note that a typical transition matrix, especially over a period 
that is not too long, is *near' / . In actual implementation of (8.1), it may give rise to 
possibly negative entries, though small in size. In this case, one makes an adjustment in the 
resulting matrix, for example by replacing such entries by 0 and decreasing the adjusted 
amount from the diagonal probability. In most cases, this procedure works well and is robust. 

We summarize various constructions made thus far. Starting with a P transition matrix on 
the risk horizon and an initial term structure of (J default probabilities we construct series of 
transition matrices for single processes on standard observation times. These are then 
translated to produce series on a prescribed set of observation times T , first for the single 
processes and then for the joint process. Of course, for single processes, the (J series thus 
obtained can be used to imply the corresponding P series if differences in their thresholds 
are provided. 



Application: We indicate how the constructed series of O transition matrices can be used in 
evaluating credit swaps. A credit swap on a corporate bord or a reference facility provides 
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proteaion to the holder of the reference bond from the ioss of value when the bond issuer 
defaults. The buyer of the contract makes a series of premium payments, contingent upon the 
survival of either the seller or the bond issuer. Here, the scheduled dates for the premium 
payments are used as observation times T and the joint process will be that of the asset 
returns of two obligors, namely the counterparty C and the issuer B of the reference facility. 

The price of a credit swap can be computed as the expected value, under (> -measure, of the 
difference between fixed rate and contingent payments. Also, when valuing the credit swap 
for risk measurement, it is imponant to consider whether the counterparty is a seller or a 
buyer, as the two situations are not exactly the opposite of each other. The payoff function for 
the contingent payments can be given generically in the form 

-^cp =Y./i-:^i:. (8.2) 

the sum is over an exhaustive set {f-:} of all possible, mutually exclusive events, is the 
indicator function of event /;, and is the cashflow, usually of what was promised in the 
contract or some fraction of it due to seller's own default, when event li occurs. Similarly, 
one has the payoff function A',,^ for the fixed payments side (buyer). 

Denoting by i^. and t^^, respectively, the random default times of the counterparty C 
(seller) and the reference party fi , then typically E is an event with a specified range for 
values of i^. and . For example, 

E = ( event that u, < < //, and r^. > //, ) 
where < //, are two specified times, usually a pair of consecutive premium payment 
times. Depending on £, is a function involving notional amount or forward values of the 
reference facility, a strike price, recovery rates of B and C , etc. and depends on whether 

SUBSTITUTE SHEET (RULE 26) 



I- 



wo 01/52121 

7j PCT/CAOO/01440 

premium payments are made in advance or in arrears. Both E and /, can be made explicit 
according to the credit swap's specification. 



The pnce of the contingem claim is the expected discoumed cashflows under the risk neutral 
pricing measure O . It is equal to 

Given a term structure ofimerest rates, the risk free discounting, c/,/, . of cashflows /, can 
be approximated reasonably in a straightforward manner. Therefore, the remaining key 
ingrediem in pricing credit swaps is to evaluate QiE), where K usually are some joim 
events peaaming to the default times x, and t„ . In practice, all the probabilit.es QiE) 
appearing in pricing formula (8.3) can be expressed in terms of the entries of the constructed 
transition matrices. We illustrate this, for example, for the following event: 

E = { event that //„ < t„ < /,,., and //„ < s ) , 
Let /., = , /e^ ) be the initial ratings of (B,C) at time ... and B, the rating of B at time 
If, (and similarly for ( \ . etc). Then, we have 

Q(E) = 0(/i,,=DX\=D\/„) 

= ic^(fV=A^.-.-AQ,=/<J/..) 

K 

where, for the last equality, we have used the Markov property for both the jomt and single 
processes. Note that the probabilities appearing in the final stage are entries of the transition 
matrices, over periods K-.^^J and [,,,.„...]. Hence QiE) are easily attained from the 
transition matrices. This completes the final p.ece of computat.on in pricing credit swaps. 
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Part 3 - Stochastic Rates 
1* Introduction 

Credit VaR model with stochastic interest and FX rates (Credit VaR II) is an extension of the 
Credit VaR model described in the Part I (referred fijrther as Credit VaR I), The chart in Figure 
6 describes the implementation of CreditVaR II model. CreditVaR II is a purely simulation 
engine and relies on generating correlated risk factors without pre-generation (like CreditVaR 
I), and full re-pricing of the portfolio for each scenario. Outputs are Value-at-Risk and 
various other statistics for portfolio as a whole as well as standalone and marginal statistics 
for each instrument in the portfolio. 



SUBSTITUTE SHEET (RULE 26) 



wo 01/5212! 



73 



PCT/CAOO/01440 



2. Simulation Methodology 



2.1. Simulation of Base Curves 



Let 0 be the current time and h be the risk horizon. Let y{t,T) be the continuously 
compounded zero coupon rate seen at time / for a bond maturing at / + T. Note that T is the 
term to maturity, not an absolute time point. Let /(/J]J\) be the forward rate seen at time 
/ for a contract starting at / + 7; and maturing at / + /;. We assume that the zero curve 
maturity buckets are the same for all currencies and all credit ratings, and we denote the 

maturity terms by « = where /..=0</, </, <...<,^. For convenience, we 

denote the zero curve as seen at time /, byy= y^) = (y(h,,,),y(h,,, ),..., y(f,,,,)) 

and the one seen at time zero by y» = X ) = W,/, XMO,/, ).-. v(0./, )) The 

zero rate for a term between two bucket points is determined by linear interpolation. (If / < /, 
or / > then the zero rate is set to .v, or >„ , respectively.) 

Our first objective is to generate samples of y, which are correlated within itself, spreads, FX 
rates and equity returns. (Note that there is a vector y for each currency. We dropped the 
subscript for currency to avoid heavy notation.) Consider the forward rates seen at h, 
denoted by x = {x„x, x,) = (/(/»./..,/,),/(/,./„/,). -./(/r, /,., ,/J) . x can be given by 



y, =>>, and X, = -^' ' ' for / >!. (i) 



Conversely, 



Written in matrix forms. 



(2) 
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(3) 



Denoting the above transition matrix by T, then x' = T-'y' where the inverse matrix is 
given explicitly as 



(4) 



The following method is proposed to generate interest rate scenarios. Let = {r",.Y",...,x^) 
= (/ (0, /„ , /, ), /(O, /(O, , . / J ) be the forward rates seen at time 0. We assume that 

the forward rates x are lognormally distributed and that the log relative changes 

^ ^^^\^^2 ^/>)= (In~p,ln^ ^n^) follow a multivariate normal distribution 

N(^,n), where =(//,, /y, ^/J are drifts and n is the covariance matrix. Therefore, the 

forward rates are guaranteed to be positive. 



To simulate the state of zero rates at time h in the ftiture (risk horizon), we first generate 
correlated normal sample ^ =(^,,4„ from the distribution as defined above, obtain : 
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= (r"exp^, x^'exp^^), then aoply equation (1) to get continuously compounded zero 

rates y . 

2.2. Simulation of Spreads 

Given a credit rating R, , / = ! - 1 , let .v, be the spread curve for over that of the 

base curve y Thus, v + .s, is the zero curve for R^ rating. The restrictions on variability of 
the spreads are the following: a) the corporate curves should increase with decreasing the 
credit rating; b) the corporate curves should not overlap for any time point / We propose the 
following scheme for generating stochastic credit spreads. In order to reduce the number of 
risk factors in the simulation and satisfy the restrictions above, we generate the spread curves 
by a combination of Monte Carlo simulation and linear interpolation. The risk factors 
involved in simulation of spreads arc: 

- Credit spread .s, (/, ) for the highest credit rating for (typically Aaa) and for maturity 
equal to the first bucket point on the zero curve (/, ) 

- incremental spreads Av (/, ) = v (/, ) - .v^ = 2 A' - I for maturity /, 

^a:-i(^) = « e., the difference between the spreads for the lowest 

credit rating for the longest and shortest maturity. 

The simulation relies on the lognormal assumption for the above risk factors. After they are 
generated as a part of a simulation scenario, we obtain spread curves by the following 
procedure: 
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1. Obtain .y,(/,) = .y, ,(/,) + ) . Lognormal assumption on Ai, (/,) ensures that spreads 
are increasing when credit rating is decreasing, i.e. .v (/, ) > v (/, ) for / = 2,.. . 

2. Obtain the .v^. JtJ = s,. x(U) ^^kJ^h) We assume that As\,J/J has lognormal 
distribution. This way we make sure that the spread curve doesn't decrease with maturity 
for any scenario. 

3. Obtain the rest of the spreads for maturity //. by partitioning the interval [o.i\..//^^] 
proponionally to .v (/,) by taking .v//^; = ^\(fy) '^^' j^^^ ^^^^ spreads that 
increase when credit rating declines for maturity th. 

4. Obtain the interior of the spread curves by linear interpolation between \, (/, ) and .v^ (/^) . 

The corporate zero curves at risk horizon are obtained by adding the appropriate spread to the 
base curves. 

2.3. Simulation of FX rates 

The FX rates at risk horizon are assumed to be lognormally distributed. 
2.4. Simulation of credit migrations 



Credit migrations and defaults are simulated similariy to the methodology of CreditVaR I 
(see Part 1). The assumption is that credit migrations are generated by the change in asset 
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value of the obligor, which is linked to industry participation and idiosyncratic component. In 
the framework of CreditVaR II, the migrations are based on the transition matrix, 
idiosyncratic obligor returns and the correlation of industry indices with the rest of the risk 
factors and within themselves. 

2.S. Statistics 

Statistics contain means and standard deviations for some risk factors (forward rates, spreads, 
FX rates) and the covariance matrix of all the risk factors including industry indices. 

3. Input Data 

Obligor data: 

Information about obligors is organized into a database containing details of their credit 
ratings, industries, and countries. 



Porlfoho data: 

Information about financial positions is organized into portfolios of exposures. It covers 
different types of instruments such as fixed income instruments, loans, commitments, letters 
of credit, etc. 

Market data and transition probabilities: 

These include yield curves, spread curves, foreign exchange rates, transition probabilities 
from one credit rating to another for different credit rating systems (Moody 8 states. Moody 
18 states, S&P 8 states) 
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Statistics: 

Statistics contain means and standard deviations for some risk factors (forward rates, spreads, 
FX rates) and the covariance matrix of all the risk factors including industry indices. 

4. Risk Statistics and Measurement 

Credit VaR II produces risk measurements that are similar to the ones produces by Credit VaR 

I. On the portfolio level, the following quantities are calculated: 

Present value and Forward value of the portfolio at risk horizon 

Distribution of the portfolio value at risk horizon and its mean, standard deviation, 

skewness and kurtosis 

Value-at-Risk numbers at various percentile levels with 95% confidence bands 
Expected Excession (or Expected Loss Given Default) of the percentile level 

On the instrument level, the following is calculated and reported: 

Present value and Forward value of the instrument at risk horizon 

Mean and Standard Deviation of the instrument value at risk horizon 

Marginal Standard Deviation and Delta-Standard Deviation (see below) of the 

instrument for a given portfolio 

Capital amount (see below) 

Optionally: Standalone Value-at-Risk numbers at various percentile levels 
Marginal Standard Deviation of instrument 
Delta-Standard Deviation of instrument 
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5. Credit Migration Model 

Credit quality migrations are modelled according to Credit VaR I methodology ( Part 1). 

6. Simulation Engine 

First, the engine prepares the necessary data for simulation. Threshold levels of standardized 
asset returns representing credit rating changes are determined for each obligor in the 
portfolio using the transition probabilities; this part depends only on the obligor s credit class 
and not on the obligor itself Relevant data is read from the database and is stored in the 
appropriate data structures in memory. 

The scenarios are generated based on the covariance matrix and other statistics of the risk 
factors that are relevant for a particular portfolio. Each scenario is composed of the base 
curve part, spread curve part, FX part and credit migration part. 

For each of these scenarios, portfolio has to be revalued. In CreditVaR I model it was done 
quickly using pre-processed prices of instrument. CreditVaR II does full revaluation of the 
portfolio, without cash flow bucketing or any portfolio compression, based on the set of 
generated risk factors. Each time a default occurs in scenario, i.e., when the sampled asset 
return value of an obligor is below the default threshold level, a random recovery rate is 
generated according to a beta-distribution whose defining parameters are governed by 
seniority associated to the obligor's instrtiments. Finally, we obtain a distribution of portfolio 
values and from it the relevant risk statistics. 
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7. Implementation 

The Global Analytics CreditVaR II program utilizes the methodology described above. The 
core computational engine is implemented in C++ Currently, interface for Windows is being 
developed, using the ADO database technology. Windows version currently utilizes a 
Microsoft Access database. ADO technology allows CreditVaR 11 to be used with the 
majority of available databases. 

It is noted that those skilled in the art will appreciate that various modifications of detail may 
be made to the preferred embodiments described throughout this specification, which 
modifications would come within the spirit and scope of the invention as defined in the 
following claims. 
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CLAIMS 

We claim; 

I. A method of creating an adjusted credit risk transition matrix of probabilities of credit 
migration comprising the steps of 

• obtaining a transition matrix under probability measure (P), 

• adjusting the (P) transition matrix to be consistent with default probabilities under 
martingale measure (0 by making a column of default probabilities in the 
transition matrix consistent with default probabilities calculated under (O), and 

• scaling the other entries in the adjusted (P) transition matrix to compensate for the 
adjustment while retaining the relative weights among non-default classes from the 
obtained (P) transition matrix in the adjusted (P) transition matrix. 

2. A credit risk model for the analysis of complex credit instrument, comprising; 
an obtained transition matrix under probability measure (P), and 
an adjusted (P) transition matrix based on the obtained (P) transition matrix and 
consistent with default probabilities under measure (0 by making a column of default 
probabilities in the adjusted (P) transition matrix consistent with default probabilities 
calculated under (0, 

wherein the other entries in the adjusted (P) transition matrix are scaled to 
compensate for the adjustment while retaining the relative weights among non-default 
classes from the obtained (P) transition matrix in the adjusted (P) transition matrix. 
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3. The credit risk model of claim 2, further comprising a Monte Carlo simulation engine to 
generate valuation scenarios based on probability distribution computed using the adjusted 
(/O transition matrix. 

4. The model of claim 2 or 3, wherein the model is implemented in computer software on a 
computer readable medium. 

5 The model of claim 4, wherein the software is for use on compatible hardware. 

6. The model of claim 2 or 3, wherein the adjusted (P) transition matrix is used to price credit 
swaps. 

7. The method of claim 1, wherein the step of adjusting the obtained (/') transition matrix is 
performed in accordance with the following: 

R = {/^, Rt: \ ' credit classes 

R,: - corresponds to default state 

' ' = - adjusted transition matrix for the time interval .[/.,,/,] 

P,.j( fnJx) ' transition probability from credit class R to credit class 




/V;(^..^) = 0 if ./ <A' and /\. 1 



^. =^V-'(/J,.,.) 
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r,= iV-'(^,) 

- obtained (/') probability of default for the credit class R, 
(from the obtained {P) transition matrix) 

- (0 probability of default for the credit class 

- a set of credit classes (can be the whole set R = {/?, R^} 
or some subset 

' ^ ' - denotes the number of elements in G ; if G is empty, put 

C=0 



G 



N() 



standard normal cumulative density function. 



8. The method of claim 1, wherein the step of adjusting the (/>) transition matrix is performed 
in accordance with the following: 

^n" = II /^"y II - adjusted {P) transition matrix for the time interval [/,„/,] for 

obligor n 

P?.j(t„,i,) . obligor's n transition probability to migrate from credit class 

R, to credit class R^ 

PK.jOo,r,) = 0 \f j<K and K'.a ('o.'.) = I 
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P"j: - obtained {P) probability of default for the obligor // over 

interval [/,„/,] if it is in the credit class at time 

H", - (C^) probability of default for the obligor u over interval 

[/,„/,] if it is in the credit class at time 

ft - instantaneous risk-free interest rate at time / 

X " obligor // stock log-return over interval ] 

- obligor // drift function of time / 

A^O - standard normal cumulative density function. 



A method of pricing a credit swap, comprising the steps of 

• Obtaining a (F) transition matrix over a risk horizon [/.„/,] and the (O) default 
probabilities over the same period 

• Creating an adjusted (P) transition matrix and (O) transition matrix over [/u,/,], 

• Building a first series of forv^ard (0 transition matrices for standard observation 
times / = from the (O) transition matrix and an initial term structure of 
default probabilities under O over [/„./,]; 
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• Translating the first series to produce a second series, having premium payment dates 
as the observation times by working with 'fractional' powers of transition matrices for 
subintervals; 

• Deriving a third series of transition matrices from the second series for co-migration 
on credit classes for pairs of obligors; 

• Multiplying cash flow for each payment date by a risk-free discount factor and by 
probability contained in the third series and summing the computed series for all 
payment dates; and 

• Subtracting the sum from a payoff fiinction for fixed payments. 

10. A method of pricing a credit swap, comprising the steps of 

• Obtaining a (/') transition matrix over a risk horizon [/„,/,] and the (O) default 
probabilities over the same period [/„./,]; 

• Creating an adjusted {/>) transition matrix and (0 transition matrix over [/,„/,]; 

• Building a first series of forward (O) transition matrices for standard observation 
times /=f,M,i^,... from the {(J) transition matrix and an initial term structure of 
default probabilities under (J over (/„,/,]; 

• Deriving a third series of transition matrices from the first series for co-migration on 
credit classes for pairs of obligors; 

• Translating the third series to produce a second series, having premium payment dates 
as the observation times by working with 'fractional' powers of transition matrices for 
subintervals; 
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• Multiplying cash flow for each payment date by a risk-free discount factor and by 
probability contained in the second series and summing the computed series for all 
payment dates; and 

• Subtracting the sum from a payoff function for fixed payments. 

1 1. A method of pricing a credit swap, comprising the steps of: 

• Obtaining a (P) transition matrix over a risk horizon [/,„/,] and the {(J) defauh 
probabilities over the same period [/„,/,]; 

• Creating an adjusted (P) transition matrix and (0 transition matrix over [/ „/, ]; 

• Building a second series of forward (O) transition matrices for premium payment 
dates from the {(J) transition matrix and an initial term structure of default 
probabilities under (J over [/„,/,]; 

• Deriving a third series of transition matrices from the second series for co-migration 
on credit classes for pairs of obligors; 

• Multiplying cash flow for each payment date by a risk-free discount factor and by 
probability contained in the third series and summing the computed series for all 
payment dates; and 

• Subtracting the sum from a payoff function for fixed payments. 

12. A system for pricing a credit swap, comprising: 

• Means for obtaining a (P) transition matrix over a risk horizon [/„,/,] and the (0 
default probabilities over the same period [/,p^]> 
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• Means for creating an adjusted (P) transition matrix and (0 transition matrix over 

• Means for. from the adjusted (P) transition matrix and an initial term structure of 
default probabilities under O over [/„,/,], building a first series of forward transition 
matrices for standard observation times / = /,./,, /j,..., 

• Means for translating the first series to produce a second series, having premium 
payment dates as the observation times by working with fractional' powers of 
transition matrices for subintervais; 

• Means for deriving a third series of transition matrices from the second series for co- 
migration on credit classes for pairs of obligors; 

• Means for discounting cash flow for each payment date over the period by a risk-free 
discount factor and probability contained in the third series, and summing the 
computed series for all payment dates; and 

• Means for subtracting the sum from a payoff function for fixed payments. 

13. A system for pricing a credit swap, comprising the steps of 

• Means for obtaining a (/') transition matrix over a risk horizon [/,„/,] and the (0 
defauh probabilities over the same period [/,„/,]; 

• Means for creating an adjusted (P) transition matrix and (0 transition matrix over 

• Means for building a first series of forward (0 transition matrices for standard 
observation times / = z,./,./,,... from the (0 transition matrix and an initial term 
structure of default probabilities under O over [/„,/,], 
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• Means for deriving a third series of transition matrices from the first series for co- 
migration on credit classes for pairs of obligors; 

• Means for translating the third series to produce a second series, having premium 
payment dates as the observation times by v^orking with TraciionaF powers of 
transition matrices for subintervals; 

• Means for multiplying cash flow for each payment date by a risk-free discount factor 
and by probability contained in the second series and summing the computed series 
for all payment dates; and 

• Means for subtracting the sum from a payoff function for fixed payments. 



14. A system for pricing a credit swap, comprising the steps of 

• Means for obtaining a (P) transition matrix over a risk horizon [/,„/,] and the (O) 
default probabilities over the same period (/„,/,]; 

• Means for creating an adjusted (P) transition matrix and (O) transition matrix over 

• Means for building a second series of forward {(J) transition matrices for premium 
payment dates from the (O) transition matrix and an initial term structure of default 
probabilities under (J over [/„./,]; 

• Means for deriving a third series of transition matrices from the second series for co- 
migration on credit classes for pairs of obligors; 

• Means for multiplying cash flow for each payment date by a risk-free discount factor 
and by probability contained in the third series and summing the computed series for 
all payment dates; and 
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• Means for subtracting the sum from a payoff function for fixed payments. 



15. Software on a computer readable medium for pricing a credit swap, comprising 
instructions to: 

• Obtain a (/') transition matrix over a risk horizon [/„./,] and the (0 default 
probabilities over the same period [/„./, ]; 

• Create an adjusted (/') transition matrix and (O) transition matrix over [/,„/,]; 

• Build a first series of forward (O) transition matrices for standard observation times 
/ =/,7,. /,.... from the (O) transition matrix and an initial term structure of default 
probabilities under (J over [/„./,]; 

• Translate the first series to produce a second series, having premium payment dates as 
the observation times by working with 'fractional' powers of transition matrices for 
subintcrvals; 

• Derive a third series of transition matrices from the second series for co-migration on 
credit classes for pairs of obligors; 

• Multiply cash fiow for each payment date by a risk-free discount factor and by 
probability contained in the third series and sum the computed series for all payment 
dates; and 

• Subtract the sum from a payoff ftmction for fixed payments. 

16. A method of modelling stochastic interest rates, credit spreads and foreign exchange rates 
comprising the steps of 

• Computing obligors- credit ratings at risk horizon oased on country/industry market 
index logreturns and idiosyncratic components; 
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• Computing base zero interest rates for all currencies and given maturities /„,/, 

from base forward rates --/(^-m/J corresponding to time intervals 

k-iJ.h 

• Computing credit spreads .v (/, ) = .s;_, (/, ) + Av, ) for maturity /, and credit classes 
/ = 2....A' - 1 (by currency), 

• Computing credit spread '\'.JfJ = ^s^,J/J-^A.s^,_/fJ for the lowest credit rating 
and longest maturity (by currency), 

• Computing the rest of the spreads for maturity by partitioning the interval 

[O.-V- /^ >>] proportionally to .v (/, ) by taking .v (f^ ) = (r^ ) ^^ -'^^^^ (by currency); 

• Computing spreads for other maturities /^., by linear or any other type of 

interpolation between a,(/,) and a',(/^) (by currency); 

• Computing the corporate zero curves at risk horizon by adding the appropriate spread 
to the base curves (by currency), and 

• Computing foreign exchange rates. 

17. A credit risk model for the analysis of portfolios of credit instruments, comprising; 

• Means for computing obligors' credit ratings at risk horizon based on 
country/industry market index logreturns and idiosyncratic components; 

• Means for computing base zero interest rates for all currencies and given maturities 

Am'i from base forward rates /(/o.O'/(V^:)- -/(^-i.^) corresponding to time 

intervals [/,,/,], - 
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• Means for computing crelit spreads -y. (/, ) = (/, ) + As, ) for maturity /, and 
credit classes / = 2...JC - l ;by currency); 

. Means for computing credit spread s,,,OJ = s,,,0,J + As,.,0,) for the lowest 
credit rating and longest maturity (by currency); 

• Means for computing the rest of the spreads for maturity /, by partitioning the 
interval [0.s^_,(/j] proportionally to s,(f,) by taking s,(lj = s,(t, )^^^^iiiid 
currency), 

. Means for computing spreads for other maturities by linear or any other 

type of interpolation between a,(/,) and (by currency). 

• Means for computing the corporate zero curves at risk horizon by adding the 
appropriate spread to the base curves (by currency); and 

• Means for computing foreign exchange rates. 



(by 



1 8. The credit risk model of claim 1 7, wherem, all risk factors (interest and foreign exchange 
rates, market index returns, idiosyncratic components) have a joint probability distribution so 
that credit migration process and interest and foreign exchange rates are correlated. 

19. The credit risk model of claim 1 8, fUnher comprising a Monte Carlo simulation engine to 
generate valuation scenarios based on Monte Carlo simulations of risk factors. 
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Obligor Data 

Credit ratings 
Market participation 
Idiosyncratic data 



Market Data 
Yield curves 
Transition probabilities 
Exchange rates, etc. 



Portfolio Data 

• Bonds 

• Loans 

• Loan commitments, etc. 




Portfolio preprocessing 

• Pre-computing exposure values 

• Compute asset correlations between the 
obligors 

• Aggregation by obligor and by seniority 




Analytic Engine 
Compute mean and standard deviation of 
exposure forward values for a given time 
horizon 

Compute marginal standard deviation and 
marginal risk by obligors 
Compute total standard deviation and total 
mean of portfolio forward value for a given 
time horizon 



Simulation Engine 

• Compute mean and standard deviation of 
exposure fonward values for a given time 
horizon 

• Compute total standard deviation and total 
mean of portfolio forward value for a given 
time horizon 

• Compute Value-at-Risk at various percentile 
levels together with confidence bands 




Final Report 

• Report the results 

• Mean and standard deviation obtained by 
two different methods should be close 



Figure 1 
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Obligor Data 
Credit ratings 
Market participation 
Idiosyncratic data 



Market Data 

• Spot base curves 

• Spot spreads 

• Spot FX rates 

• Transition probabilities 
' Market data statistics 

(mean + covariance matrix 
of base curves, spreads, FX 
rates and industry indices) 



Portfolio Data 
Bonds, Loans 
Floating Rate Notes 
Interest Rate Swaps 
Loan commitments 




Simulation Engine 

* Simulate scenarios consisting of correlated credit 
rating migrations (and defaults), corporate zero 
curves at time horizon. FX rates at time horizon 

• Re-price portfolio for each generated scenario 
Obtain the distribution of portfolio value at time 
horizon 

Obtain Mean, Standard Deviation, Value-at-Risk 
(with confidence bands). Expected Excession of 
portfolio value at time horizon for given percentile 
levels 

Obtain standalone and marginal statistics for each 
instrument 



^ Y 

Final Report 
• Customized reports of the simulation results 



Figure 6 
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